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ANALYSIS OF RADIATIVE HEAT TRAITSFER FOR 


LARGE OBJECTS AT METEORIC SPEEDS 
By Kenneth K, Yoshikawa 
Ames Research Center 


SUMMARY 


An explicit closed-form solution has been obtained for radiative heat 
transfer to a body in flight in the earth's atmosphere at speeds such that 
radiation is the dominant mode of heat transfer. The solution was attained 
by assuming a gray-gas radiator and then linearizing the radiative transfer 
equation. Approximations to the one-dimensional flow equations are developed 
which allow the solution to be applied to the stagnation region of a shock 
layer on an axisymmetric body. Because of its characteristics^ the linearized 
solution is readily applicable to multilayer heat transfer encomtered during 
very high speed flights^ when the interactions among the cold air ahead of 
the shock wave^ the heated air behind the shock wave^ and the layer of 
ablation gases adjacent to the body surface must be considered. 

The development of the first-order linearized solution for radiative 
heat transfer is described^ and some applications are presented and discussed. 
The results demonstrate the effects, on radiative and convective heat transfer 
to the body surface, of shock-layer energy loss by radiation, absorption of 
radiation within the multiple layers, blockage of radiative heating by blowing 
of ablative vapors, and the absorption and emission of radiation by the body. 


INTRODUCTION 


For atmospheric entry at meteor velocities, radiation is the dominant 
mode of energy transfer. Many investigators have emphasized the importance of 
understanding and predicting the radiant heat transfer to bodies at these 
speeds (e.g., refs, 1-3) • These investigators have pointed out that energy 
loss by radiation from the shock layer and the reabsorption of radiation in 
the shock layer itself, in the free stream, and in ablation products will 
affect both the radiative and convective heat transfer to the body. 

Both simplifying flow assumptions and intricate numerical techniques have 
been used to solve the gasdynamic problem. Examples of the flow assumptions 
that have been used are: one-dimensional porous flow (ref. k) and a linear 

velocity distribution and constant pressure in the shock layer (ref. 5) ■ 
Reference 6 showed that the mass flow distribution on the stagnation stream- 
line, rather than the velocity distribution, is relatively insensitive to 
radiative heat transfer. Examples of mmierical techniques are local 
similarity (ref. 7 ) and the integral method (ref. 8). 



Most investigators have emphasized the effects of radiation cooling 
rather than self-absorption on shock-layer structure. The inclusion of self- 
absorption into the problem introduces additional complexities that require 
fiarther simplification of either the equations or the forms of their solu- 
tions. These s amplifications generally involve the substitution of an expo- 
nential fimction for the exponential integral that appears in the solutions of 
the transfer equation, and/or the iise of successive approximations employing 
combinations of simple fimctions (refs. 9*“13)« 

Substituting the exponential function is the simplest means of evaluating 
reentry radiation, including self -absorption, with reasonable accxaracy. 

At relatively high free-stream density, where the radiation mean free 
path is comparable to or smaller than the body size, one can assume that radi- 
ation flux is absorbed by the media just ahead of the shock front (preheating). 
The structure of the strong shock front has been investigated (refs. l4-17) 
for plane shock flow of a perfect gas. Such analyses explicitly consider 
thermal radiation effects on the structure of the shock wave rather than the 
influence of preheating energy on radiative heating. 

The ablative mass loss from large meteors entering the atmosphere was 
considered (ref. l8) assuming a porous-flow analogy and the Rosseland approxi- 
mation. The increase in total shock stand-off distance due to a sizable 
amount of ablation or injection has been shown, both theoretically and experi- 
mentally, to obey a linear correlation formula (refs. I 9 and 20); this pro- 
vides a basis for calculating radiative interaction between shock and ablation 
layers . 

Our primary interest here is to obtain the simplest form of an analytical 
solution for radiative transfer in the shock layer, especially explicit 
expressions for radiative net fluxes at the body and at the shock. This solu- 
tion includes the effect of self -absorption in a strong shock layer, the 
interaction of thermal radiation with the shock layer, the preheating zone, 
and the ablated vapor layer. The results will be used to relate the 
convective heating rates in adiabatic and nonadiabatic flow. 
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radiative flux absorbed by preheating zone (^gp - 
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direction cosine or viscosity 

quantity defined by equation (F3b) 

order of Bessel function (eq. (35b)), or frequency 

parameter, mnax^ 

density 

absorption coefficient, cm“^ 

Stefan-Boltzmann constant 

optical thickness defined by equation (2a) 
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optical thickness at interface 
optical thickness at wall 

radiative blockage function at wall (eq. (107)) 

r 

radiative leakage fimction, — — 
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thermal function defined in eqmtions (l8c) 

Exiler's constant (0-57721 . • . ) + 2 Zn 2 (eqs. (55)) 
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shock wave due to preheating effect 
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ANALYSIS 


The eqiiation of radiative transfer will be combined with, a simplified 
gasdynamic energy equation for one~dimensional flow^ and the thermodyiamic 
properties of air will be introduced in such a way as to make the resulting 
equation linear. Solutions for this linear eqixation will be obtained for 
plane shock flow and stagnation stream-tube flows. The extension of the tech- 
nique to multiple gas layers (e.g.^ a zone of preheating or of injected or 
ablated species) will be considered. The calculation of radiative heat- 
transfer rates at boundaries will be emphasized. However^ auxiliary results^ 
such as the distribution of enthalpy and temperature in the shock layer and 
the reduction of shock-layer thickness as a result of radiative transfer will 
be obtained. Details of the analyses are deferred to appendixes wherever 
possible . The notation used follows closely that of reference 5 or of 
reference 21 with minor differences* 


Assumptions and Geometry 

The following assumptions are basic to the analysis: 

1. One-dimensional flow (plane shock flow and stagnation stream flow) 

2 . Local thermodynamic equilibrium 

3 . Gray-gas approximation 

k. Viscosity and heat transfer by thermal conduction are neglected. (How- 
ever^ the effect of radiation on convective heating will be estimated in 
a later section.) 

5* In the solution of the transfer equation, the exponential integrals are 
replaced by exponential functions. 

6. The term T^ is assumed to be locally a linear function of enthalpy, and 
a thermal function is defined which incorporates this relation (locally 
linear thermal relation) . 

Figures l(a) and l(b) show the geometry used in the analysis. 


Basic Radiative Equation 

If the absorption coefficient of the gas is assumed to be independent of 
frequency (gray gas) and the gas to be in local thermodynamic equilibrium, the 
radiative transfer equation for plane parallel flow is 

^ P + I = B (1) 
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where I is the specific intensity (rate of radiative energy flow per unit 
area and solid angle), and p is the direction cosine relative to the free 
stream (upstream direction is positive) . The optical thickness t and the 
Stefan-Boltzmann function B are defined by 


-I: 


r = I pK dy 
o 


(2a) 


B 


-/ 
w o 


dv = -r T^ 


(2b) 


where By is Planck’s fmction and k, the mass absorption coefficient. 
The formal solution to equation (l) is 


= l“(0)e 


-(V-t^) 


^ r - 


-(T-t)/-p 


dt/-|i 


(3a) 


+ + -(tw-t)/p r'^w -(t-x)/p 

I (T,h) = I (T^)e + J 


Be 


dt/h ( 3 b) 


where I is the specific intensity in the upstream direction (positive p) 
and I~ in the downstream direction (negative p)j is the optical 

thickness at the downstream boundary. 

The rate of radiative energy transport per unit area (fig. 1(b)) is 


q (t) = 2jt 
q'*’(T) = 2jt 

It follows that 

q"(T) = 2qgE3(r) + 
q*^(T) = 2q(^3(T^ - 


f I (t,p)p dp 
^ o 

(4a) 

I'*’(t^p)p dp 

(4b) 

2it f BE2(t - t)dt 
^ o 

(5a) 

t) + 2it ^ BEa(t - r)dt 

(5b) 
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•where q“ = rtl“(0), ) , and the exponential integral function E_ 

(refs. 9 and 22) is ^ 


= r ® ^ dp. = f e dx 

(6) 

^0 

The net flux transferred per unit area q is 

9.(t) = q'*'(T) - q (t) 

(7) 

If the dimensionless notations 

r = F(t) = q(T)/2rtBg 

(8a) 

r; ^ <i;/2«Es = 

(8b) 

II 

CQ 

OJ 

III 

(8c) 

T = T/Tg 

(8d) 


are introduced into equation (7)^ we obtain 


r(T) = -2R^E3 (t) + 2R^3(t^ - t) 




T^EaCt - T)dt 


(9a) 


The subscript s refers to position immediately "behind shock wave. The first 
and second temis in equation (9st) represent the radiative fluxes at t from 
the upstream and downstream boundaries^ respectively^ attenuated by the gas 
between t and the boundaries . The third and fo-urth terms represent the 
radiative fluxes at t from the gas between t and the upstream and down- 
stream boundaries^ respectively^ attenuated by self -absorption. At the 
upstream and downstream boundaries^ equation (99') reduces to 

Tw 

fg = r(0) = -2RgE3(0) + 2R^3 (t^) + f T%(t)dt = T+ - Rg (9b) 

^ O 


= r(T^) = - 2 RgE 3 (Tj + 2R^3(0) - l^E 2 (T^-t)dt = Rv- (9c) 

where E 3 (o) = l/2. 
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Eqmtion ( 9 a) is applicable to mixLtiple radiating and absorbing layers 
as well as single layers (e-g.^ a preheating zone in front of a shock wave or 
a layer of injected or ablated gases from a s-urface) . However, the exact cal- 
culation of fluxes is complicated by the interaction between layers. To 
simplify the calculation an approximation is introduced whereby the actual 
boundary flioxes are replaced by "equivalent wall” fluxes . This approximation 
concentrates the effects of other layers at the boundaries of the layer for 
which the solution is desired and allows equations (9) to be applied to each 
layer of the multilayer separately. The details of the procedure are 
presented in appendix A. 


Pseudo One-Dimensional Flow Model and Thennal Function 

In hypersonic flow, as the gas travels downstream from the shock wave, 
its enthalpy decreases as a result of radiative loss to the surroundings . 

Since this study is concerned with flight conditions for which radiative heat- 
ing is the dominant mode of energy transfer (e.g., for a meteoric object 
entering the atmosphere), the energy equation in the shock layer may be 
written 


pv dh = dq 


(10) 


where heat conduction and changes in the gas kinetic energy have been 
neglected. The mass flow along the stagnation streamline is expressed in 
terms of optical thickness. First, we introduce the mass flow distribution 

X = pv/pboVoo (11a) 

Many authors (e.g., refs. 13 , 23 ) have assumed a linear relation between 
mass flow distribution and physical thickness 

X = (1 - y/L)^ (lib) 

where 7 is a constant (zero or one) that characterizes the flow. Flow with 
7=0 corresponds to plane shock flow such as in shock tubes, piston problems, 
and porous type flows; flow with 7 = 1 corresponds to stagnation streamline 
flow. For present purposes, it is convenient to replace physical thickness 
in equation (lib) with optical thickness. Thus, we assume in this paper that 

X = (1 - i-J (11=) 

where 7 has the same meaning as before . The validity of this replacement 
is discussed in a later section, where it is shown that equation (lie) is as 
good an approximation to the actual mass flow distribution as equation (lib) 
or better. With the mass flow distribution of equation (lie) the energy 
equation becomes 
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PcoVoo(i - <aii = 


(12) 


The thermodynamic variables are now expressed by a dimensionless 
function (thermal fmction) F so that 


dF = 


PoqVoo ^ 

2aTg^ 


(13) 


or 

- F = (hg - h) = ^ (1 - h) ilk) 

‘^crXs 

where is the ratio of the black-body radiation at the shock wave to the 

total flow energy. 


^a ~ Poo^oo^S 

h = h/hg 


After integrating this last equation for a constant pressinre process (or 
constant density),, one can readily obtain (ref. 4) 


where 




N 

n=o 



1 


(15a) 


ai = -8aT„^/p V c-Q ) 

S' 00 00 Ps ( 


(151) 


etc • J 


The variations of h and F with temperature are shown in figirre 2 for two 
reference conditions. 


It can be seen from figure 2 that the temperatiure variation can be 
further simplified for wide ranges of enthalpy as 

= 1 - a(Fg - F) (15c) 

if a is taken from the average slope of (3 t^/8f) at the designated 
reference point. ThuSj 
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a 




(I5d) 


Appendix B gives a more detailed 

From equations ( 13 ) and (12) the 
related to the radiative function by 

dF = — 

0 


discussion of the thermal function, 
dimensionless thermal function is thus 



Linearization of Radiative Transfer Equation 
The differential form of eqmtions ( 9 ) is 




2R^E2 \ 


(t„ - t) - 


2T^ 


T^Eidt - t| )dt 

o 


( 17 ) 


Substituting equations (l5c) and (l6) into eqxaation (I 7 ) } we obtain a linear 
integrodifferential equation with a singular kernel 



+ 2acp - a 



cfEi(Ie - e'l = f(®) 


(18a) 


with an initial value of 


Cp(Tw) = 0 


(l8b) 


where 

cp = Fg - F 

e = Tv - T > 

f = (1 - 2Rg)E2(T„ - e) + (1 - 2Rv)E2(e) 


(18c) 


No attempt is made to find a general solution for equation (l8a) except for 
7=0, for which several numerical calculations have been carried out by a 
successive iteration method merely to compare it with the approximate 
solution . 
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Approximate Solutions for a Single Layer 


In the present paper an approximate solution to the transfer equation is 
given explicitly in analytic and closed form when the e 2 q>onential integral 
function E 3 (t) is replaced hy an exponential function- We introduce the 
commonly used approximation 

Es(t) ~ me (19) 

where choice of m and n will depend on the optical thickness at the wall 
(t^) of the plane layer ^ and will he constant for given t^. The values of 
m and n are discussed later; meanwhile^ the solution is carried out without 
assigning values to these constants. 

It follows from equation ( 19 ) and the properties of the exponential 
integrals that 


= -E 3 (t) - mne 


-nr 


( 20 ) 


The integral equation (9a) then becomes 
r = -2mRse”^’^ + 2mRwe 


-riT , _ _-n(Tw-T) _ ^ r-T n(T t) 

J o 


dt + m J''"' dt 

( 21 ) 


After differentiating the above equation twice ^ we obtain the following 
second-order nonlinear differential equation: 

d^r ^ ^ df^ ar o 
+ 2mn - n r = 0 

dT^ dT 


( 22 ) 


(Equation (22) is equivalent to the Milne-Eddington approximation except with 
different constants.) From eqmtion (l6) the derivative of the temperature 
becomes 


dT^ = 

dT 


S(Fg - F) S(Fs - F) X dT 


(23) 


With the simplified equations (lie) and (l5c) and the transformation 


0 = Tw - T 


(24a) 


d0 


dr 


(24b) 



A linearized form of equation (22) is now obtained as 


^ ^ - n 


d 02 '\9y d 0 

with bouadary values^ from equation. ( 2 l)^ 


T = 0 


P W 

fg = - 2 mRs + ^ J ^ 


’’’w _4 -nt 


e dt 


(25a) 


(2513) 


= “-2inRge 


-nr. 


w 


2 mRy - mn / 
^ o 


-4 -n(T^t) 


T e 


dt 


(For consistency^ E 3 (o) = m is used instead of l/2.) 

The general solution of equation (25a) has the form 


( 25 c) 


r = ciX^(ne) + C 2 Y^(ne) ( 26 a) 

where ci and c^ are constants and X and Y are functions to be determined. 
In eqmtion (26a) X and Y for 7 = 0 are the ordinary exponential functions^ 
and for 7 = 1^ the modified Bessel fimctions. Solutions will be obtained in 
explicit and closed form for these two cases. 

It will be shown that the solutions of equation (25a) lead to the fol- 
lowing important linear^ relations between the boundary values Rg and R^ 
and the fluxes at the boundaries: 

fg = Mq + MiRg + MaRw (26b) 

= Ms + MiRg + MsRw (26c) 

Expressions for the M functions will be given in terms of d and optical 
thickness 

Plane shock flow (7 = O) The differential equation (25a) for 7=0 
becomes 


d^r 

dT^ 


+ 2mnd — - n^P = 0 
dT 


( 27 ) 


4liis is a consequence of choosing a linear thermal function (eq. (l5c)) 
and does not require the exponential approximation (eq. ( 19 ))* 
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Its solution is 


r 


-nPiT -nP2(''’w“''') 

c le + Ca© 


where 


Pi = njx + 




+ 1 


Pa = -mx + + 1 


Flxuces at the boundaries are 


r 


s 


~nPaT-w 

Cl + ca© 


= Cie 


-nPiT, 


w 


+ ca 


To evaluate the constants ci and ca^ combine equations (29) with 
equations (25b) and (25c) 


-nPaT^ 

:i + Ca© = -2mRg + 2mR^e + mn / T e dt 

o 


-nPiT^ -nr. 

Cl© + Ca = -2mRs© 


^ + 2mR^ - mn r ^ T^e 

J n 


Tw ^ -n(Tw-t) 


dt 


where (t) may be expressed in terms of Ci and Ca by substituting 
tions (28) into equation (16) and integrating between 0 and t. Thus 
equation (l5c). 


Fg - F = i (1 - T^) = ci(l - e + Ca© 

a 


Substituting eqtiation (31) into equations (30) gives two simultaneous 
algebraic equations for ci and ca of the form 


KiCi + KaCa — K5 


K3C1 + KiCa - Ks 


where 


(28a) 

(28b) 

(28c) 

(29a) 

(29b) 

(30a) 

(30b) 

equa- 

with 

(31) 

(32a) 

(32b) 
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Ki = l+ma*^l--e 




.-nTw _ 1 

1 + Pj 


1 - e 


-n(x+Pi)T^ 


]} 




-nP2*^W/ 


1^1 + loa 'll 


Kg = e 7 1 + ma *^1 - e 


1 


1 - P2 L 


1 - e 


-n( 


}) 


K3 . m [l - 


K 4 = 1 - ma 


r -nPaT^^j^ _ ^ 1 - e \ 

1 1 + P2 L J i 


Ks = 2m ('i - Rg) - 2me - 


H.) 


Ks 


= -2ni (1 - Rw) + 2me"“'^'l(i - Rg) 


J 


Cl = 


ca = 


(-KaKe + K4K5 )/(KiK 4 - K^Kg) 
-(-KiKe + K3K5 )/(KiK 4 - KaKs) 


( 32 c) 


( 32 d) 


With some algebraic manipulation, equations ( 29 ) and (32) can be ejspressed in 
the linear form: 


Ts = Mq + MiRg + MaR^ 


( 26 b) 


^ + M4Rg + MsR^ 


( 26 c) 


where the M values depend only on m, n, a., and and are given by 
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M - i m Cl - 

^ A 


(Ke + K4) - + K3) 




Mi = - 2 inj 


-UT^/ -nPaT-rr 

e IK2 - e ^Ki J - ( K4 - e 


-nPaTw 


K3 


Ma = ^ 2mJ 


Ka - Yk 4 - 


Ms 




A 


;l 


(Ki + Ks) - e ^(Ks + K4) 


\ (33) 


M4 = ^ 2 n^e - e - (ks - 


-nPiT„ 
e ^Ki 


Ms = ^ 2mj 


-nPiT. 


Ki - e ■ ‘ ’^Ka ) - e "( Ks - e "' -"■ ”K4 




A = K1K4 - KaKs 


J 


The terms Mq and M3 are fluxes from the shock layer without boundary con- 
ditions; Ml ani^ M4 show the effects of reflection or emission from the shock 
■wave (at t = O); similarly, Ma and Ms show the effects of reflection or 
emission at t = For the present case, 7 = 0 , it can be shown that 

Ml = -Ms • 

Stagnation streamline flow (7 = l) •- The differential equation ( 25 a) for 
7=1 becomes 

^ _ 2mnaTy ^ _ n^F = 0 (^4) 

d02 e d0 ^ 

which can be transformed to the modified Bessel equation of order v (refs. 2 h 
and 25); its general solution is 

r = e^[cilv(ne) + CaKv(n0)] (35a) 

where and are modified Bessel f'unctions of the first and second kind, 
respectively, and 

V = innaT^ + l /2 ( 35 b) 
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The constants ci and are evaluated by a procedure similar to that for 
plane shock flow. (Appendixes C and D give detailed evaluations of properties 
of the Bessel f\mctions required for the present analysis.) 

Since dF = df/(0/T^) y for fy = r(0)^ it follows that 


de 


nTwf 


v-i 


(36a) 


where 


Tv-i ^[cil^_^(n0) - CaK^,_j^(n0)] (36b) 

By integration one obtains 



nT 


w 


CiC^_j0) 


C2C_(^_^)(0) 



Thus, 


Fs - F = nx^ - C2C_^^_^)(rJ - j^ciC^,_j0) 

where 


C2C_(^_^)(0) 



(37a) 






(e) 



t I (nt)dt 
v-x 


t’' \^_^(nt)dt 


(37b) 


(37c) 


Then K takes the following forms: 
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where 


where 


“IlT^ T n 

Ki “■ “HinciT^ ^v-i 2 2 ^w^v-i 




-riT 


w^ 


1 * 


Ka - mnaT^e + _ Z_^ - - + - e Z_^ 


1 1 

Kg = -innaT^C^_^ + - e + - e 


Ki = 


Ks = 


1 I -riT^ , 

^ 2 ® 2-v - 2 ^ ■^wZ_(v-x) ^ 2 ^-v 


2m ('I - Rs) - 2me - R^ 


Ke = -2m (^1 - + 2me - Rg 




(38a) 


“^v-l “ 

C/ \=C/ 

-(v-l) -(v-l)'' W' 

Z_v = Z_^(t^) = T^^K^(nT^) 


z"",, ^ Z_„(0) = i f^) r(v) 


-'__v ' 2 


(38b) 


j 


r(v) in eq-uations (38^) a gamma function. Bomdary fluxes are 


fg - CiZ + c^Z 


-V 


= "az!, 


(39a) 

( 39 b) 
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The relation between fliixes and boundary values is again 


Tg = Mq + MiRg + MeR^ 


= Ma + M<tRs + MsR„ 


where M now has the following forms: 


Mq = I m (1 - e [(Ka + K4)Z^ - (Ki + K3)Z_v] 




Ml = 4 2d 


Mp = — 2mJ 
A 


-nT. 


w,. 


-nT 


w 


e "Ka - K4 - e Kj. - K3 Z_ 


— nT Tj \ / — nT -y 

Ka - e K4 ) Z^ - (Ki - e K3 j Z.^ 


Ms = m fl - e 


^ (Ki + Ks)Z*^ 


M4 = 


2m ( e - Ks ) Z*, 


1 / w \ * 

Ms = i 2m ( Ki - e K3 ] Z_^ 


A = K1K4 - K2K3 


J 


Ttie same physical interpretation can he ascribed to these values of M 
those in equations (33) • From equation (35a) it follows that the flux 
gradient at the wall is 


^w = 0 


The enthalpy gradient at the wall from equations (36) depends on v 


(26b) 

(26c) 


(i^O) 


as to 

(41) 
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for V < 1 (tea) 

for V > 1 (teb) 

•where 

dO 

The physical interpretation of this dependence of enthalpy gradient on 
V -will be considered in a later section. 

Approximate Solutions for Optically Thin and Thick Layers 

In this section^ the solutions to equations (26) are simplified for 
optically thin and optically thick shock layers. Only the major results of 
these simplifications are presented here. Details are presented in 
appendixes . 

Plane shock flow (7 = O) The following are the results for two 
limiting cases. 

Optically thin layers (t^ « l): If the exponentials are approximated 

by Taylor series^ the quantities K in equations (32c) become^ to order 

T ^ 

Ki =: 1 + i mnaPiT^s 

K 2 — 1 - ^ 

Ks ^ 1 - nPiT^ + I n^Pi(Pi 

Ki - 1 + ^ 

K 5 - 2 m(Rs - R^) + inn(l - 

Kg -2m(Rg - R^) - imi(l - 

A £1 n(3i + P2 )t^( 1 - nPaT 

With the relation derived from equations (28) 
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I nil nil 



Pi - Pa = 2ma 


ikk) 


The M values then directly follow from equations (33) and are 


Mq - mnTT,) 


1 - ^ (1 + 2ma)T^ 


Ml -=: -2m + 0 (t^^) 


Me 2m(l - nT^) 


Ms - -mnx^ 


1 - i n(l + 2ma)T^ 




M 4 — - 2 m(l - nx^) 


Ms - 2 m + 0 (x^ 2 ) 




(i^5) 


The fluxes at the "boundaries for m = l/2 and n = 2 are 


Tg — [1 — (1 + a)x^]x^ — Rg + (1 — 2 Xv■)R•^J^ 


(46a) 


^ -[1 - (1 + S)t^]tw - (1 - 2 x^)Rg + R^ (46h) 

Here^ we have omitted second order terms from the boundary fluxes because^ in 
practice^ Rg and R^ will usually be small. The choices of m and n will be 
justified later. 


Optically thick layers (t^ » l): For large optical thickness^ K takes 

the following forms (the second and third identities for Ki follow from 
eq.s. ( 28 )): 

Ki ^ 1 + mPi/(l + Pi) = K 4 - (1/2) (1 + Pi) 


K 2 - ( 1 / 2 ) (1 - fe)e ~ 0 

K3 — -md. 

K 4 — 1 + nia.[l/ (1 + 3a) ] 

Ks - m(l - 2Rg) 

Ke - -m(l - 2R^^) 


(47) 
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The determinant A becomes 


A (1 + 

Values of M follow directly 
Mq - 2m/(l + Pi) 

Ml - -km/il + Pi) 

Ma — [4m/(l + Pi)][(Pi + Pa)/(l Pi)]® 
Ms - - 2 mP 2 /(l + Pi) 

M 4 — -4m(l - Pa)/(l + Pi) 

Ms - ^m/(l + Pi) 

Net fluxes at both boundaries from equations (26) are 



( 48 ) 


(49) 


(50a) 


r 2m^ + 2(Pi - l)Rg - 2piR^] (50b) 

^ 1 + Pi 

Flux variation and temperature variation for 1 « t « are 

r = (51) 

1 + Pi 

= p^ + (1 - P2)e“''‘^'-'' + 2(1 - Pa)(l - e“''^^'')Rs (52) 

The gas temperature at the wall (t = » l) is 

V “ + ^•^2(1 - Pa)Rs + 2(1 - Pa)Rv (53) 

Effective surface temperatures that emit equivalent radiant fluxes from both 
boundaries are 

(Tef)s = (2rs)'-^^Ts (54a) 

(Tef)w = (2rw)^/^Ts (54b) 

2^1- 





It is noted from equations ( 50 ^) and ( 53 ) that radiant flux toward the wall 
for an optically thick layer with no boundary effect is always greater than 
the local black-body radiation at the edge of the layer. This was also shown 
in references 4 and 12. 

Stagnation s treamline flow (7 = l) Appendixes E and F give details of 
the derivations for « 1 and » i^ respectively, since the manipula- 
tion of the Bessel functions becomes complicated. The following are the 
results . 


Optically thin layers (t^ « l): For v = [(I/2) + mnar-^j-] fixed and the 

argument of the Bessel functions approaching zero, the values of K in 
equations (38a) are 


Ki 

K2 


K3 

K4 


Ks 


Ks 

A 


nr^^[l - mnaT^jQi 

[1 - (2mndcpQ + nT^^^°)T^]Q2 
Z 




n 


+ 2mnaJ t. 


w 




-n(mna)T ® 


(1 - 2mnacPoT^)Q2 


Z* 

-V 


-2m(Rs - Rv) + mn(l - 2R-w)'T-y^ 
- 2 m(Rg - Rv) - mn(l - 2 Rs)tv 
nT^2(i - 2mndcPo'’'v)QiQ2 j 


( 55 a) 


where 


^o 

^o 

Qi 


Q2 


2 Zn 2 + 0.57721 (where second term is Euler's constant)! 



\ ( 55 b) 


j 
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r(v) and r(l - v) in eqmtions (55i>) are gamma functions, 
are 


Mq mnT^ 1 - | n(l + 


Ml ~ -2m[l - 0 (t^2)] 




The vaJLues of M 


Mg 

Ms 


:r 2m(l - nr^) 


rr -naiTU 


1 - 


1 

2 



n(l + 4mx)T^ 


(56) 


M<t if -2m(l - nT^) 

Mg ~ 2m[l - 0 (t^2)J j 

The boundary fluxes for m = l/2 and n = 2 are 

Tg 2f '^y[l ~ (l 2a)T^] - Rg + (l - 2 t^)r^ (578’) 

- (1 - 2t^)Rs + % (5Tb) 

The thermal function is 

F, - F = (l/y{l - [1 - - nvoT/*'° - [1 - v„(n + a?o)T/^ojR^ 

- [1 - Vo(n - 2 'Po)t/'’°]R„} ( 53 ^) 


Which, to order can be f\irther simplified as 

Fg - F - (l/d){l - [1 - (T/T^)f'°}(l - Rs - Rw) (58b) 

From equation (l4) , and letting Rg = R^ = 0 in equation (58b) the enthalpy 
profile within the layer becomes 

h = 1 - (Aa/a){l - [1 - (t/t^)]^^°] (59) 

which is similar to the expression in equation (15) of reference 23 for 

« 1 . 
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The temperature distribution follows directly from equation (l5c) 

T [1 - (60) 

which is equivalent to that of reference 5 for « 1. 

It is interesting that the "boundary fluxes for plane shock flow 
(eqs. (46)) and stagnation flow (eqs. (57)) differ only in second order of 
Ty. for Ty. « 1 . 

Optically thick layers (t-^ » l); Asymptotic expressions for K and M 

are: 

Ki I (1 + Pi)Zv 

- I (1 - 

K3 - I (1 + Pi)(l - Pa)Zv 

K4 i (1 

Ks ~ 2m (^ - Rg^ 

Ks - -2m - Rv^ 

A - ^ (l + 3x)(l + -J Jtv)Z-^;Z*y 
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2m 




o 1 + P] 


Ml 


- 4 m 

1 + P; 


Mb - 


4 m Pi + Pa ^-v 
1 + Pi 1 + Zf V 


Ms - -2m 


_J2 

1 + .JUv 


M4, :r - 4 m 


1 - Pa 

1 


Ms - 


4 m 

1 + sjHv 


J 


(62) 


Thus, the "boundary flioxes are 

Tg - ^ (1 - 2Rg) (63a) 

1 + Pi 


- -2m _ [1 + 2(Pi - l)Rg - 2PiR^] (63"b) 

1 + ^TtV 

As noted from equations (62) and ( 49 ), fluxes with no "boundary effects (Mq 
and Ms) are nearly independent of the constant n (Ms in eqs. (62) is very 
weakly dependent on n) . The major contributing factor to the flux calcula- 
tion comes from the constant m. Note that the asymptotic values of for 

stagnation flow decrease inversely as the square root of whereas for 

plane shock flow reaches finite asymptotic values . 

Asymptotic radiant flux variations for 1 « t « are 

r ^ Mq( 1 - 2Rs)Zv(0)/Zv 

(64) 
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and the thermal f met ion becomes 

Fg - F ^ Mod - 2Rg) 


1 - ^ 


"'’W Zy(0) 


9 Z 


V j 


. r,(i - 


Finally^ the temperature distribution is vrritten 
1 - Mq (1 - 2 Rg)aj 


_ Zy(0 ) 

0 Z 


V J 


( 65 ) 


= Pa + (1 - + 2(l - p2)(l - e-‘^Pi'^)Rg 


( 66 ) 


The asymptotic values of the radiant fmetions above yield identical resTilts 
to those for 7 = 0 which shows that the asymptotic values of radiant 
properties are independent of velocity profile for 1 « t « 


is 


2 


For an optically thick layer ^ the Rosseland approximation of net flux 


r ^ 



dr 


Applying this equation to equation (66)^ we obtain 

r = - I npi5Mo(l - 2Rs)e"’^^^^ 


Which differs from eq-uation ( 64 ) when combined with equation (63a) by a factor 
( 2 / 3 )nPid. If the Rosseland approximation is to apply^ higher derivatives of 
temperature distribution must vanish^ or nPi « 1 . Since d is a finite 
constant^ it is required that (2/3)(n3i)d « 1. Now it is obvious that the 
coefficient (nPi)d cannot be mity^ in fact^ Pi > 1; consequently^ the 
Rosseland approximation for a high enthalpy flow field with radiation^ in 
general^ is not applicable. This also implies that agreement with Rosseland 
approximation for an optically thick layer in a flow field is not a property 
of the asymptotic solution (this approximation can be commonly used for non- 
flow cases) and it may sometimes lead to incorrect solutions if proper care 
is not exercised. 
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Calculation of Boundary Conditions 


Boundary values Rg and Ry may be expressed in terms of emissivity and 
reflectivity of boundaries since the boundary values introduced in 
eq-uations (8b) and (c) can be written® 


‘Is = 






+ (external fluxes from upstream) 


t flux from ”1 flux reflected 

shock wavej from shock wavej 


(6?a) 


w 


= 




fliuc from 

1 

1 

wall 



flux reflected 
from wall 


+ (external fluxes from the wall) 


(6Tb) 


In dimensionless form, equations (6ja) and (b) may be written (for no external 
fliuces) 


R = — T'^+rF'*' 
•^s 2 s -^s s 


“4 m 

Rw = + r„r 


2 w 


W -w 


(68a) 

(68b) 


where Tg = Tq/^s “ '^w/^s * equations (9) and (26) it follows 

that 


F'*’ = Mf. + (1 + Mi)R„ + M2R. 


(69a) 

(69b) 


F^ = -Ms - M4Rg + (1 - Ms)R;^ 

Combining equations (68a) and (b) with equations (69a) and (b) gives 

i {[1 - r„(2m - M5 ,)]€sT/ + raU2€„T„'‘} + ro{ [1 - r„(2m - Ms) - r„Ma%} 




{[1 - rs(2m + Mi)][l - r^(2m - Ms)] + rgr-^^M2M4} 


(70a) 


^Radiative flux will be reflected back into the shock layer, for example, 
just prior to an impact . Radiative flux from a surface placed Just ahead of 
the shock can be treated similarly. 
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^ {[1 - rg(2m + Mi)]e^T^^ - rJM^egTs^} - v^[[l - rg(2m + Mi) jMa + rgMoM^} 

Rv = 

{[1 - rs(2m + Mi)][l - r^(2m - Ms)] + rgr^MaMi} 


( 70 b) 

Bomdary temperatures Tq and may be taken as effective temperatures so 
that they could include fluxes impinging from external soTjrces added 
independently from the system. 


Selection of Constants m and n 

Replacing the exponential integral function by an exponential function is 
essential to the present analysis and is commonly used as an analytic 
approximation for radiative transfer problems. Here the relation 

E 3 (t) = me ( 19 ) 


is introduced where m and n are constants chosen in some plausible fashion. 

Although many different values for m and n are found in the literature^ 
depending on the applications^ the first constant^ m^ is not arbitrary; it 
must be l /2 to satisfy the asymptotic solutions for both optically thin and 
thick layers. This is clear for optically thin layers since 

E3(0) = ni = I (71) 


For optically thick layers^ m must also be l/2^ as will be shown in Results 
and Discussion where exact and approximate solutions of the equations are 
compared. The exponential constant n is less restricted^ and proposed 
values of n between 3/2 and 2 are found in the literatiire. The value 3/2 
is associated with Eddington ^s approximation and has been preferred for astro- 
physical problems (ref. 2l) . In this paper, values of n which depend on 
total optical thickness are \ised to provide correct results for both 

optically thin and thick layers. The criterion is chosen to match the area 
under E 3 with the area under the exponential function: 


(72f 


"^Quick convergence results if equation ( 72 ) is rearranged as 
G = - 1 + 2 n[(l/ 3 ) - E 4 (t^)] and solved for e = 0 . 


f ^ E3(T)dT = m f 


w -nT 
e dT 


or 


m(l - 

n = 

I - E4(t^) 
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Thus, for Tv -»• 00 , E 4 (tv) 0, and n = 3/2 "when m = l/2. For Tv -*• 0, 
-Es(0) = £ 2 ( 0 ) = mn and n = 2 when m = l/2- A set of exponential constants 
is tabulated in table I for various values of and for m = l/2. 

This completes the solution of the radiative transfer equation in terms 
of optical thickness - 


Conversion of Optical Thickness to Physical Thickness 

To facilitate application of the solution to physical problems the 
optical thickness must be converted to physical thickness. Physical distance 
into the shock layer is related to optical depth by 

dy = dr/pK (73a) 


Integrating eq.uatlon (73a) between 0 and gives the nonadiabatic standoff 
distance L as 


L = 



dr/pK 


( 73 ^ 3 ) 


For stagnation flow^ Goulard (ref. 5) has shown that the nonadiabatic distance 
into the shock layer ^ y^ is related to the adiabatic distance Y by 

dY = (p/pg)^^^ dy = (p/pg)^^^ dr/pK (74a) 

Therefore^ the adiabatic standoff distance is 

Lo = <1 t/pk (7i^-b) 

o 


From equations (73^') s*3id (7^3*) *the ratio of nonadiabatic to adiabatic standoff 
distance is 


L/Lo = <i(Y/Lo) 

o 


(75a) 


An approximate relation between adiabatic distance and optical thickness 
(ref. 5) is 


d(Y/L^) d(r/T^) = dT 


(75b) 
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Incorporating this approximation into equation ( 75 a) gives the following 
simple expression for standoff distance^ which is consistent with the present 
approximation for mass flow^ eqioation (lie). 




^o 



( 75 c) 


In order to use eqmtions ( 7 ^ 1 ^) and ( 75 c) the density ratio (p/pg) and 
absorption coefficient (pK) can be expressed in terms of power series of 
enthalpy ratio (ref. 26), and the shock-layer thickness readily evaluated by 
numerical calculation. Finally, the Hayes-Prob stein formula relates Lq to 
body radiiis and density ratio across the normal shock wave (ref. 27) by 



(76a) 


where 


k - 



and an average value is 



for — ~ l6 
P 

oo 


(76b) 


(T6c) 


(See also ref. 28.) 


Approximate Solutions for Multilayers 

Absorption of radiation by the free-stream (preheating) If appreciable 
amounts of radiation are absorbed by the free-stream media just ahead of the 
strong shock front, there is an important effect on the structure of the 
shock layer and the radiative heating to the body is changed significantly. 
(Emission effect from upstream is somewhat similar to increasing the reflec- 
tivity of the shock front.) Following references l 6 and 29 the present paper 
assumes negligible change in kinetic energy and nearly constant density ahead 
of the shock wave^ (as opposed to nearly constant pressure in the shock 
layer); the thermal function for air in the preheating zone is therefore 
determined by a constant density process. Furthermore, the flow in front of 
the shock wave is planar and the analysis previoiisly presented for 7 = 0 can 
be applied to the preheating region. Therefore, with the "equivalent wall" 
approximation (appendix A) for the boundary conditions at the shock wave, the 
flow can be separated into two regions with variable boundary conditions: 


5 

These assumptions are considered in more detail in reference 29. 
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Region * is the preheated free stream with P “ Poo (constant) y = 0. 

Region I is the shock layer with p Pg (constant) and 7 = 1 . 

In the preheating zone, the energy eqixation is 


Pa.Voo(ll* - hoo) - (l* - Qoo) ^ 0 ( 77 ) 

.n the preheating zone and qoo is the 
(sketch (a)). Physically, this model 
is adopted because, in real air, only 
the far ultraviolet portion of the 
flux leaving the shock wave will be 
absorbed in the free stream near the 
shock layer. Consequently, in the 
present analysis the portion of the 
radiation leaving the shock front 
(q^) is allowed to escape from the 
system. No attempt is made to com- 
pute q^, but it simply becomes a 
parameter in the analysis. 

In the shock layer the energy 
equation may be written from 
equation (lO) as 

pv d(hsp - h) = d(qgp - q) (78) 

where qgp is the net radiative heating at the shock wave due to preheating. 
Equation (78) shows that preheating changes the enthalpy level in the shock 
layer. This increase is proportional to (qgp - since energy balance at 
the shock wave (7=0) is 

PocXoo (tloo 2 ^ ~ “ Poo^oo (^^sp 2 ^S ^ ~ * 3 .sp 


where q* is the net radiant heat flux i 
radiant heat flux escaping from the zone 


I ^co ~ 0 

r 


Vco 




Rf 




Rs 


sp/ 


Too T* 

Preheating Shock wave 
edge 



Wall 


Sketch, (a) Preheating and shock layer. 


which can be rewritten as 


PooVoo(l^Sp - llg) = <lsp - < 3 oo 


(79a) 


since 


^ ^ - i V 2 

» 2 00 2 s 


Thus, the increases in enthalpy immediately ahead and behind the shock are 
obtained. 


Prom equations (78) and ( 79 a) three types of energy profiles can be 
established . 
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(a) qgp q«, 

(^) <3oo ^ ° 

(c) <lsp > <loo > 0 

Case (a) is for no preheating and is the case previously discussed. Case (h) 
is for infinite layers that reach adiabatic equilibrium Rankine-Hugoniot 
temperature as ^ oo (for 1 « t « t^) . Case (c) is a general one, and 
contains strong absorption, or a preheating effect and will be discussed here. 
In dimensionless form equation (793') becomes 

Fsp Tgp - (79b) 

where fgp = Q.sp/^^®s* From equation (l^c), the corresponding texnpera,tvore 
rise behind the shock wave is 



+ d(r 


sp 


Foo) 


( 80 ) 


Temperature in the preheating zone is usmlly relatively small compared 
to that in the shock layer; thus emission from this zone may be neglected. 

(An extension of this part of the analysis which inclijdes emission from the 
preheating zone is given in appendix G.) With this assumption and f;jrther 
assuming negligible wall emission (Rg = R^ = O) , we have, for given and a, 
from eq-uation ( 26 b) 

^sp = ^pMo (8la) 

—4 

where Tgp must be inserted to account for the higher temperatinre behind the 
shock wave due to preheating. Similarly, from equation ( 26 c), the flux at the 
wall is 

= TgpMs (8lb) 


We now define 


(jj = 



(82) 


Combining this definition with eqmtions (8o) and (8l) gives 

n _ Mo 


sp 


wp 


1 - oMo^ 


Ms 
1 - oMo 


10 


(83a) 

( 83 b) 
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(83c) 


(1 - <^)Mq 
1 - aJ^w 


Since emission in the preheating zone has been neglected, the transfer eq.m- 
tion (9a) in this zone takes on the particularly sin 5 )le form 

r* = rgp2E3(T*) (8i+) 

We use equation (84) to define so that = t* when f* = Too- Hence 

J’co = rsp2E3(Tj (85) 

and is fixed by w through the relation 

Es(tJ = I (1 - u)) (86 ) 

Finally, the enthalpy distribution in the preheating zone (from eq. (7?)) is 

— ~ ^ = Aa(r* - rj (87a) 


which at the shock front becomes 

= >^(rsp - r„) (87b) 

“s 

where h^ is enthalpy immediately ahead of the shock- Pressure and tempera- 
ture in this zone are now obtained directly from the thermodynamic chart for 
constant density (p* = p^) . The distribution of thermal properties in the 
shock layer with preheating may be found by adjusting the thermal function 
previously found (eqs. (l4) and (l5c)) to account for the higher ten 5 )erature 
due to preheating immediately behind the shock wave. If the normalization 
temperature for the thermal function F and the f met ion f is changed from 
Tg to Tgp, then the previous analysis of the shock layer, equations (35) and 
( 37 ) # nmy be applied directly. Therefore we define a new thermal_function in 
the shock layer, F'^", normalized with respect to T^p by Ft = F/T^p, or 

(Fsp - F)'^ = (Fgp - F)/Tgp (88a) 

This new thermal fwiction is related to the original thermal function Fg - F 
(without preheating) by 

Fg ~ F = (F 3 - Fgp) + (^sp 

= -(Fsp - 0 + Tjp(Fgp - F)t (88b) 
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The enthalpy and temperature distributions in the shock layer obtained in a 
similar manner are 


h = 1 + AaUFgp - rj - T^p(Fsp - F)^] (89) 

^ - ^(Fsp - F)^] (90) 


Injection and ablation layers .- It is assi^med that a secondary layer of 
vaporization materials from the wall is separated from the primary shock layer 
hy an interface. Experimental observations support this ass-umption (refs. 20 
and 30 ) • Figure 1 indicates the geometry with a stagnation point in the flow 
at the interface. The analysis of the vaporization layer is similar to the 
analysis for a shock layer in that constant pressure with 7 1 (stagnation 

flow) is assumed. As in the shock layer ^ it is assumed that the thermal prop- 
erties in the vaporization layer can be correlated by a linear equation relat- 
ing and enthalpy. The reference temperature for this correlation is 

taken as the temperature of the vapor at the wall. It is further assumed that 
the injected gases are relatively cool compared to the air in the shock layer; 
therefore^ the emission from the injected vapors may be neglected. An exten- 
sion of the analysis in which this last assumption is relaxed is presented in 
appendix H. 

We define a radiative blockage function as 

^wb _ ^wb 
^wo r^o 

where is radiative flux penetrating to the wall through the vapor layer 

and fvo is the radiative flux for the same flight conditions without a vapor 
layer but, for consistency, with the same wall temperature with and without 
vaporization . 

The optical thickness of the vapor layer is measured from the wall to the 
interface and conditions in the vapor layer are denoted by the subscript b . 
Thus, = 0 at the wall and t^-|^ is the optical thickness of the vapor 
layer at the interface . 

If the emission from the vapor layer is assumed to be small compared to 
the emission from the shock layer, it follows that the properties in the shock 
layer are not changed significantly by the presence of the vapor layer. 
Therefore, from equations ( 9 )^ 

Fwb - (92) 

Consequently, 

( 93 ) 


( 91 ) 
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The physical thickness of the vapor layer depends on the mass flow rate 
of the vapors; in adiabatic flow the thickness of the vapor layer is 
correlated by a simple expression (refs. 20 and 3l) • 


h. ~ 

Lo ■ 


1 


+ k-b 





where = Lq + Lbo-’ total distance between shock wave and wall; 

kb - 1 for a spherical nose (ref. 20); and are the molecular weights 
of air and vapor layers; Lq and are the thicknesses of air and vapor 

layers for adiabatic flow; pb is the density of the injected gas at the wall; 
and fv is the blowing ratio ( Pb''''b/Poo^oo) • 

The correlation law (eq.. (9^)) further provides a simple relation between 
vapor layer thickness and the blowing parameter; thus. 


^o ^ /Ps Ho 
Lo ^7 Pco Mb 




(95a) 


^o 



(95b) 


where kg - 0.75 for spherical noses (see eqs. (?6)) and R is the radius of 
the spherical body. The nonadiabatic thickness of the vaporization layer is 
approximated by equations (T^b) and (75c) in a manner similar to the 
calculation of nonadiabatic shock-layer distance. 


For forced injection of vapors from the s\rrface, the blowing rate param- 
eter f^ is given a priori and the solution is then straightforward. For 
ablation, however, the blowing rate parameter f^ depends on the heat reach- 
ing the wall F-wb# 8-C-<3- the solution is therefore coupled. We ass\mie that the 
mass loss rate of ablation vapors per unit area is given by 


dm ‘^'wb 

- 37 ^ Pb^b = — 


(96) 


where ^ is the heat capacity of ablation products. The net radiant heat 
through the entire system at the wall is 

l^b ” ^‘^s^^wb ~ '^Vb'lwo (97) 


thus, the blowing parameter, f-w = 


Pb'V’b/PocVoo^ becomes 


,4 
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( 98 ) 





^wb 


where \ = (2oTs‘^/pboVootis) as in equation (l^) . The coefficient a-^ in the 
thennal function for the ablation layer is coupled with the blowing parameter 
f^^ since 1^-5 eqimtion (98) depends on co^ for given optical thickness 
and boundary conditions as given below: 


where 


and 



“ PooVoJlw 

^ a(T/Ty)^ 
\Bh ^ (h/hv) 


( 99 ) 


With equations ( 95 )^ (98)^ an^ ( 99 ) ^^^.d the solutions previously pre- 
sented for 7 = 1, including equation (T^b), an iterative procedure is used 
to determine f and f^ and the ablation vapor-layer thickness Lt)o (see 
appendix H) . ^ 

A simpler and useful approximation for may be obtained if it is 

further assumed that the absorption coefficient (P|^^-b) the ablation layer 
is constant and that the blowing parameter f^ is very small. 

Then, by the definition of optical depth, 

^wb ^^b^ (100) 


where L-^, is nonadiabatic vapor layer thickness, equation (lOO) may be 
rewritten as 


"^vb . 2 ^o-^voS B 


(101a) 


where 
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With equation { 93 ) > % becomes 

_ -i^orwog^B 

'Vb ® 


( 102 ) 


From equations (9?) (98) and by tbe definition of B in equations (lOlb), 

can be written 


\b 



(103) 


Equations (102) and (IO3) may be combined to eliminate B and therefore to 
obtain O^b (°^ vice versa) for ablation. For B « 1 (small ablation), 
equation (102) may be expanded in the first order of power series and combined 
with equation (IO3) to obtain the approximate formula 


®wb 


1 

1 + ^Qf^Qg 


(104) 


where n has been set equal to 2 and Lb/L-^Q = 1. Equation (104) is useful 
in that the absorption coefficients of the shock layer and vapor layer (which 
are generally unknown) are combined into a ratio . This allows parametric 
studies for assessing the effects of radiative blockage . 


RESULTS FOR A SINGLE LAYER 


Calculations of radiative flux at the boundaries and of profiles of 
radiative flux, enthalpy, and temperature in the shock layer have been made 
for selected flight conditions . The equilibriim air properties used in making 
these calculations are given in appendix I. The res\iLts are used to assess 
the effect of radiative energy loss on shock standoff distance, the effect of 
preheating and absorption in ablation layers on the radiative heat transfer 
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to the body, and the effect of radiative energy loss on the convective heat 
transfer. Where possible, the effects of the assumptions on the present 
results are considered by making comparisons with other solutions . 

Comparisons With More Exact Analyses 

Plane shock wave flow (7 = O) To evalxmte the error introduced in the 
present solutions by replacing the exponential integrals with exponential 
functions, equation (l8a) (the integrodifferential equation) was solved 
numerically^ and this solution compared with solutions of equation (25a) (the 
differential equation obtained by using the exponential approximation) . 

Various commonly used values of the constants m and n were used in these 
calculations. The results are shown in figure 3(0.) for 7 = 0, where and 
Fg, the radiative heat to the wall and free stream, respectively, are plotted 
as functions of t^. This figure shows that the exact and approximate solu- 
tions agree well for m = l/2 and n varying from 3/2 to 2 (but constant for 
a given value of t^) as provided by equation (72) • As mentioned in the 
analysis, m is not arbitrary but should be l/2 as required by the asymptotic 
values for optically thin and thick layers and by the equivalent wall approx- 
imation. This becomes particularly apparent (fig. 3 (si)) Fg as optical 

thickness is increased since the results for values of m different from 
1/2 deviate from the correct asymptotic fliuc. The other constant n is less 
restricted, except that it should take on the value of 2 as approaches 

zero because it must take mn = 1 for ^ 0 . Note that the approximate solu- 
tion for m = 1/2 and n = 3/2 shows rather good agreement with the exact solu- 
tion except, as expected, for small optical thickness. 

Stagnation flow (7 = l) A comparison similar to that of figiire 3(a) 
but for 7=1 has not been made primarily because of numerical difficulties 
associated with the singularities in equation (l8a) . However, Howe and 
Viegas* (ref. 7) elaborate viscous solutions are available for comparison 
with the present solution. They do not assume a linear thermal function, a 
linear mass flow relation, or the exponential approximation (as does the 
present work) . However, their solutions are for a gray gas, with the same 
absorption coefficients as used in the present work. Furthermore, the solu- 
tions of reference 7 si^re for relatively small values of and thus energy 

loss by radiation is emphasized rather than energy redistribution by self- 
absorption. The comparison is presented in figure 3 (b) which shows very good 
agreement between the two treatments except, of course, near the wall where 
viscoias effects predominate. The approximate boundary-layer edge, 6, is 
indicated in figure 3(b) as estimated from the empirical equation (ref. 32) as 

3 ( 105 ) 


sin the iteration process for equation (l8a) the approximate solution is 
substituted as an initial input function. This reduces the difficulty of con- 
vergence associated with the equation and is very effective in controlling 
oscillations of the solution. 

kl 
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The stagnation streamline mass-flow distribution uised in the present 
work is shown in figure 3(c). The calculated mass-flow distributions for 
vario\is optical thicknesses are plotted as a function of physical depth 

into the shock layer and are compared with both a linear eq.mtion and the 
more realistic quadratic equation (nonadiabatic compressible flow) presented 
in reference 6. This figure shows that the present results agree well with 
the quadratic eqviation for most optical thicknesses but approach the linear 
equation for small optical thickness. As a first approximation, the present 
mass-flow distribution (eq. (lie)) provides the advantage of a closed form, 
simple solution for the radiative heat transfer. 


Radiative Flux Results 

The principal results of this report are illustrated in figures 4 and 5. 
The radiative fluxes at the shock wave Fg, and at the wall, as determined 
from the present closed- form solution are plotted as f met ions of in 

figmes 4(a) and 4(b) for a temperature Tg = 15,000° K and pressures p^ = 1 
and 10 atm, respectively. Results are presented for both plane shock flow 
and stagnation streamline flow; the corresponding flight conditions are shown 
on the figure. As mentioned in the analysis, radiative fluxes for plane 
shock flow and stagnation flow differ only to second order in optical thick- 
ness when is small. Furthermore, the asyn 5 )totic values of Fg for 

large Ty are the same for both flows. The asymptotic absolute values of 
^w for stagnation flow, however, decrease inversely as the square root of 
Ty, whereas Fy for plane shock flow decreases very slowly to its finite 
symptot ic value . 

The results shown in figures 4(a) and 4(b) for stagnation flow are 
replotted in figure 4(c) where the ratio of nonadiabatic to adiabatic radia- 
tive heat transfer Fy/Fo is plotted as a fmetion of the ratio of adiabatic 
radiative heat transfer to flow energy (Aq) . Existing numerical machine solu- 
tions are also shown on figure 4(c) . For Pg = 1, where absorption effects 
are small, the present method agrees well with the results of references 7 and 
8 except when Aq is small. This can be anticipated because for small Aq 
and fixed free-stream conditions Fq must be small, which in tvim corresponds 
to small nose radius. Therefore, the Reynolds numbers for these conditions 
are low and the shock layer is primarily viscous. Thus, the temperatme pro- 
file deviates sharply from that for a radiating field alone, and the present 
method overestimates radiative heating. It is probable that the present 
resvilts for pg = 1 in figure 4(c) represent a correlation curve for all 
flight conditions as long as Ty is small and thus absorption effects are 
small. (This has also been pointed out in ref. 8.) However, when self- 
absorption becomes pronomced, as is the case for Pg =10, Fy/FQ deviates 
from this correlation c\n’ve. Thus, the correlation curve should not be vised 
for large values of 7\q, since large values of Aq eventually correspond to 
large Ty. 

Results similar to those presented in figure 4, but for a much higher 
velocity (30 km/sec) and for stagnation flow only, are presented in figure 5 . 
The ratio of radiative heat to total flow energy q/(l/2) is plotted in 

figures 5(a) and 5(l>) for Pg = 1 and 10 atm, respectively. The same trends 



are evident in these figures as in figures 4 (a) and 4 (b) . Figure 5 (c) shows 
parameters similar to those in fig\jre 4 (c) except that the parameter is modi- 
fied by A* = (vj,gf/Voo)Ao^ where the reference is chosen for the flight con- 
dition, Pg = 1 atm and = 15 km/sec. An e3q)lanation of the correlation 

is that the modified parameter A* depends weakly on the flight velocity 
since 

Ps'^s “ Ps ^ ^s^ and p^Voodg “ Psn/^ 

or 

- ^00 


Thiis, parameter Ag = Wref/^oo^^Q becomes relatively Independent of the 
flight speed. For the higher velocity flight condition, and for both pres- 
s-ures, values of A* less than about 1 correspond to small values of 
Thus, the results correlate well with each other and with the correlation 
curve from figure 4 (c) . 


Shock-Layer Profiles 

Enthalpy distributions in the shock layer for various values of Ty are 
shown in figures 6(a) and 6(b). For large Ty the major change in enthalpy 
occurs relatively near the shock wave; for small the major change in 

enthalpy occurs near the wall. Enthalpy, when plotted in terms of optical 
thickness t (rather than f = T/xy), is always higher for larger Ty. 

Because of the nearly constant enthalpy near the wall for large Ty, 
radiation heat transfer to the wall for this case (l.e., for large bodies or 
high free-stream densities) is nearly black-body radiation where 

is the local temperature near the wall . Such radiation is proportional to 
local enthalpy, as can be seen from the thermal function, and therefore is 
proportional to a reduced velocity vj defined by 


h 




(I06a) 


Ihe absorption effect thus reduces radiative heat by an amount which 
corresponds to a velocity change similar for a conical body with apex angle 


= 



(I06b) 


An important parameter in the analysis for 7 = 1 is v (the order of the 
modified Bessel function where v = (I/2) + mnaTy) since the character of the 
solutions depends on 'whether v is greater or less than 1 . In figure 7 ^ 



profiles of T, r, and dr/dT are presented for v > 1 and v < 1. The 
temperature distribution for v < 1 shows a very sharp drop near the wall and, 
in fact, has infinite slope at the wall, whereas the temperature distribution 
for V > 1 has a finite slope near the wall (eq.s. (i|-2)). The temperature 
itself at the wall is greater for v > 1 than for v < 1. The flux distri- 
bution shows no infinite slope; in fact , the derivative of the flux is zero at 
the wall regardless of the value of v (eq. (4l)). 

In figure 8, the variation of the gas temperature at the wall, Tg (with 
optical thickness, t^), shows considerable differences between stagnation flow 
and plane shock flow. These differences occur because for stagnation flow the 
flow time (and thus the emission time) near the stagnation point is infinite. 
For stagnation flow of an optically thin layer (t^ « l) the gas does not 
absorb strongly enough to block the radiant energy, and the edge temperature, 
therefore, is relatively low. As the optical thickness increases, as a result 
of increasing either the absorption coefficient or the body radius, Tg 
increases with increasing since now the gas near the wall absorbs flux 

from upstream more effectively. A fijrther increase of t^, however, cavises a 
gradual drop in Te, and Tg finally vanishes as approaches infinity. 

This occurs because the gas particles then have sufficient distance to travel, 
th\is time in which to cool, even in a strongly absorbing layer. Because the 
velocity approaches zero at the wall in stagnation flow, there is a coupling 
effect between self-ab sorption and emission time for the gas near the stagna- 
tion point. Plane shock flow does not show this coupling effect (as illus- 
trated in the same figure) becavise the velocity behind the shock is 
approximately constant and therefore the emission time is very short compared 
to that for stagnation flow. In actxial flight, the boundary layer plays a 
significant role near the wall, so the trend of Tg described above exists, 
but is interrupted by heat conduction in the boundary layer before the Tg 
reaches its final value by radiation alone . 

The effects of absorption in radiating flow discussed above will also 
occ-ur for the case of nongray gas radiation where cooler gas near the stagna- 
tion point (with or without boundary layer) is heated by the absorption of 
U-V radiation from the rest of the shock layer. In other words, the optical 
thickness is in effect increased locally, and thus will act to prevent a 
sharp drop of edge temperature . 


Shock Standoff Distance 

Figure 9 presents the ratio of nonadiabatic to adiabatic shock standoff 
distance; this ratio illustrates the effect of density increase due to tem- 
perature drop by radiation loss . The results of reference 7 are included in 
this figure for comparison. The curves of L/Lq flatten out as the body 
radius increases and absorption effects become significant. For small bodies, 
the shock standoff distance depends primarily on radiation loss and is less 
sensitive to the effects of absorption. 


kk 


APPLICATIONS TO MULTILAYERS 


Blockage of Radiation by Injected or Ablated Vapors 

The effect of injected or ablated vapors on the heat that reaches the 
*wall from an intensely radiating shock layer will be considered; convective 
heat will be assimed negligible compared to radiative heat. The radiation 
blockage parajneter (defined in eq. ( 91 )) is the ratio of the net radia- 
tion flux reaching the wall with blockage to that reaching the wall 

without blockage ^wo • It is assumed that the intensity of the shock-layer 
radiation is so high compared to that of the ablation layer that the major 
effect is absorption of shock-layer radiation and therefore emission from the 
ablation vapor is negligible. The emission term for the injected vapors 
becomes comparable to the absorption term when for the vapors becomes 

very large; however^ that occurs only when values of a*3re small and thus 

radiative heating at the wall has already been greatly reduced, (it should 
be noted that emission from the ablated vapor increases the net radiative 
heating at the wall.) Furthermore, convective heating can increase the 
temperature and consequently the emission from the ablation layer. 

Figure 10(a) presents the radiative blockage fmction for injec- 

tion, which is similar to the blockage function for convective heat. The 
solid line on the figure was obtained by considering only absorption in the 
injected vapor, equation (l02) , while the dashed line was obtained by con- 
sidering both absorption and emission with various optical thicknesses of the 
shock layer, t I. A rather good correlation curve resulted for these calcu- 
lations. It can be seen that the effect of emission is negligible for 
'^wb calculations were obtained with carbon as the injected vapor. 

Fig-ure 10(b) presents the radiative blockage function, for ablation 

as a function of the ratio of absorption coefficient in the ablation layer to 
that of the air behind the normal shock wave for body radii of 1 and 
10 meters. Increasing the body radius reduces the radiative blockage param- 
eter for a given absorption coefficient. To block radiative heat effectively 
a small body (l m) requires higher absorption coefficients than a large body 
because the small body receives less radiative heat; consequently it has a 
relatively thin ablation layer that will not block as much radiative heat. 
Radiative heat is sharply reduced as the absorption coefficient of the abla- 
tion vapor increases because the optical thickness in the ablation layer is 
increased by two factors: 

(1) Higher absorption coefficients 

(2) Increased physical thickness due to the absorption effect (i.e., 
the ablation layer is heated by absorption, and accordingly is 
thicker) . 

In connection with the flow field, temperature near the interface 
changes sharply; but it is continuous for stagnation type flow (7 = l) since 
net flux and velocity at the interface are continuous. The theoretical 
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interface ten^jeratuire for a sim il ar sitmtion in plane shock flows (7 = O) is 
discontinuous because the velocity profile is discontinuovis . 


Preheating Zone 

As the flight velocity approaches meteoric speed, most of the radiation 
from the strong shock layer is in the U-V region and can be strongly absorbed 
by free-stream air ahead of the shock. The absorption of the ultraviolet 
creates a preheating zone which raises the temperattire both ahead of and 
behind the shock wave, thus increasing the radiative heat to the wall. To 
represent the actual (wavelength selective) absorption by the free stream 
with a gray-gas analysis, a portion (Poo) of the flux entering the preheating 
zone from the shock layer (fgp) is allowed to escape from the system (as 
described in the analysis). The escaping fraction is then ^>oo = ^ot/^sp 'which 
can be considered a radiative leakage function at the edge of the preheating 
zone. The fraction absorbed by the preheating zone is 1 - ^>oo. The Increase 
in heating to the wall is given by ^yp where ^>vjp = ^’wp/^wo^ ratio of 
radiative heating to the wall with and without preheating. 

Figure 11 is a plot of $wp versus body radius R, with 4bo as a 
parameter. (Emission from the preheating zone was neglected in figure 11.) 

For $00 = 1 there is no preheating effect and = 1. For ^00 = 0 all 

radiative flux is trapped in the preheating zone, and there is a sizable 
increase in radiative heat with increasing body radius. This increase is 
pr ima rily due to the increase of temperature behind the shock wave. For a 
given (escaping radiation fraction), the larger bodies undergo a greater 
Increase in radiative heat at the wall because the preheating zone can absorb 
more from the larger radiation associated with the thicker shock layers of the 
large bodies . Figure 11 also suggests that in actual flight a body can 
receive more than half the energy radiated from the system (Qoo "*■ *lwp) • How- 
ever, as the temperature ahead of the shock becomes con^arable to that behind 
the shock, emission and other neglected effects (i.e., kinetic energy changes 
and heat conduction) may become important. 

The approximate combined effects of both preheating and ablation on heat 
transfer to the body can be given as 

^ ^b (107 ) 

preheating effect "1 fablation vapor effect 
on nonablating bodyj |_ without preheating 

if absorption is the major mode of energy transfer in the ablation layers . 


Effect of Honadiabatic Flow on Convective Heating 

While it is beyond the scope of this paper to analyze in detail the 
changes in convective heating that occur as a result of nonadiabatic radiative 
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flow^ it is possible to make a simple estimate of the effect by modifying 
somewhat the criterion (first introduced by Goulard^ ref, 33 ) of taking the 
enthalpy at the edge of the viscous boundary layer as the driving enthalpy for 
convective heating. (Later ^ Thomas (ref. I3) suggested taking the enthalpy at 
the edge of the gas layer ^ h^ (the enthalpy corresponding to T^) , as the 
driving enthalpy.) It is reasonable to assume that convection becomes impor- 
tant in energy transfer at the boundary-layer edge and finally becomes com- 
pletely dominant in changing the fluid enthalpy at a certain distance into the 
boundary layer. The assumed distance at which convection becomes dominant^ is 
6/2^ where 6 is the boundary-layer thickness given by equation ( 105 )- How- 
ever^ the additional loss of enthalpy by radiation in the transition layer 
should be considered where radiation and conduction are both important. To 
approximately account for this additional loss of enthalpy by radiation we 
will use the enthalpy at 6/2 computed for radiation alone ^ ^6/2^ driv- 

ing enthalpy by convective heating. Thus^ the convective heating qc 
reduced by 


^^5/a 

^o 


(108) 


where q^^^ is convective heat without radiation effect. Figure 12 shows the 
ratio Q.c/^co ^ function of t^. The results agree well with the results 
from reference 7*® 


CONCLUDING REMAKKS 


The gray-gas approximation was used and the radiant heat-transfer equa- 
tion was linearized to obtain analytical closed-form solutions for the radia- 
tive heating of a body for optical thicknesses from zero to infinity. It was 
shown that effects of absorption and energy loss are important. By virtue of 
the linearized solution, radiative heat transfer in multiple layers with 
various boundary conditions is simplified considerably. Some important char- 
acteristics of nonadiabatic flow were computed simply from adiabatic rela- 
tions, for example, shock standoff distance in terms of density ratio and 
convective heating in terms of enthalpy ratio. A radiative blockage function 
for gas injection, similar to the convective blockage function, was introduced 
and computed for high-speed flows over ablating bodies in which radiant heat 
predominates over convective heat. It was also shown that the preheating 
effect in high-speed flight may contribute substantially to the heat transfer 
to a body. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif., 9 ^ 035 , March 20 , I967 

129-01-08-11-00-21 


"^The precise value was discussed in reference 29- 

®The asymptotic value of Iq/Icq approaches zero slowly since Tg 
approaches zero . 
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APPENDIX A 


RADIATIVE TRANSPORT EQUATION IN DOUBLE OR MULTILAYERS 


A preheating layer with a shock layer and a shock layer with ablation or 
vapor injection from the body sinrface separated by an interface represent a 
typical two-gas layer problem. The radiative equation within the shock layer 
for this case becomes (from eq. {9o.)) , for t < Tj^ 


r(x) = -2RgE3(T) + 2R^3 (Tw ~ ~ J '*'S ^ T"^E2(t - r)dt 


(9a) 


which, in terms of each layer thickness, becomes 


r(T) = -2R E3 (t) + 2R^3(t^ - - t) - / T^E2(t - t)dt 


r.Tx 


+ / ^ T^E2(t - T)dt + / ^ T^E2(t - r)dt 


(Al) 


For T > , 

r(r) = 2RgE3(T - Tj_ + T^) + 2 R^ 3 (t^ - t) 


^ T^E2(t - t)dt - 


T^E2(t - t)dt + 


'w 

T^E2(t 


- x)dt 


(A2) 


Equivalent Wall Approximation 

It is usef\£L here to introduce the concept of an equivalent wall approxi- 
mation, that is, the assumption of a thin wall that intercepts all incoming 
fluxes on this boundary, and reemits the same number of photons in the same 
direction across the wall. Thus, mathematically, equations (Al) and (A2) 
change form slightly and become, for t < t^. 


r(r) = -2RgE3(r) + 2 RiE3(t^ - 



T^E2(t 



T^E2(t - r)dt 


(A3) 
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for T > 


- 

r(x) = -2R2E3 (t -T j^) +2R^E3 (t^-t) -J T^E2(r-t)dt + J T^E2(t-T)dt 


where 


r ^ - 

Ri = 2R^3 (Tw - Ti) + / T^E2(t - T^)dt 


(a4) 

(A5) 


R2 = 2RsE3(Tj^) + / T^E2(xi - t)dt 
J o 


(a6) 


Qioantities Rx and represent equivalent boundary conditions which include 
the effect of radiative flux from other gas slabs on the assixmed thin wall. 
Consequently^ if equivalent wall conditions are assumed^ radiative transport 
of double layers is separable and can be considered as two independent single- 
layer problems with variable bomdary conditions. The physical interpreta- 
tion given by equations (A3) to (a 6) can be justified mathematically as well. 
By Taylor’s expansion^ for ^ 


EsCt) = E3(t^ + t - t^) 

= E3(tj^) - (t - T^)E2(Ti) + . . . 

and 

E3(t^)2E3(t - T^) = E3(tj^) - (t - Tj^)2E3(Tj^) + . . . 

- E3(t^) - (t - Ti)E2(Tj^) + . . . 

for the moderate range of (exact at t = Tj^) . Since the major contribu- 
tion of the exponential integral function E 3 (t) comes when t and the 

magnitude of the difference of the above two functions becomes small for 
larger t 


1E3(t) - E3(Tj^)2E3(t - Ti)| « E3 (t^) 

for the entire range of t. Thus^ the above substitution of E 3 (t) by 
E 3 (Tj^) 2 E 3 (T - T^) will provide a good approximation for all values of Tj^. 

In fact^ when the integral exponential function is replaced by the exponential 
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function of constant m and n^ the original transport equations (9) and 
modified equations (A3) and {Ah) become identical; thus^ from substitution of 


E3 (t) = me 


-nT 


and its derivative (Ea(T) = -EsCt)) 


E2 (t) = rone 


-nT 


it follows for T > Tj^ that 


„ , s -n(T-Ti) -nTi 

EsCt) = me ' e 


E3(t^)2Es(t - T^) = 2m^e ^e 


2 -nTj_ -n(T-Tj_) 


Thus^ 


E3(t) = E3(Tj_)2E3(T - Tj_) 


(A7) 


when m = l/2. Similarly, for t <Tj_, 


E3(Tw - t) = E3(t^ - Tj_)2E3(Tj_ - t) 


(a8) 


With the exponential approximation substituted in the integral part of 
equation (Al) 


-4 , , r^^-4 -n(t - t) 

T E2(t - T)dt = / T mne 


dt (t < Tj_) 


and 


r w -4 

2Es(Tj_ - t) / T E 2(t - T^)dt = 2me 

T-; 


-n(Ti-T) p'^W ^ ) 

/ T mne ^ dt 


. 2m r " ' '">dt 


(A9) 
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One can readily see that these integrals are identical if m = l/ 2 . As 
discussed previously in the approximate solution, the constant m is found 
to he 1/2 from other criteria. 


Equivalent Boundary Values 

Since the transport equation in double or miiltiple gas layers can be 
divided into two or several separate layers, corresponding boundary conditions 
on the equivalent wall can be treated, from eqimtions (26), if regions I and 
II are considered separately. Thus, region II can be considered entirely 
separate from region I so that its flow field and notation become consistent 
with that of the first layer (see sketch (b)). The superscripts I and II 
denote properties in the respective regions. 


Rs R| |R2 

Rw 

< 


rI 

Rw 

r nr 











iw 

1 

1 



y®() 



1 

0 Tj 

1 

V 

and c 

) 

T 

(Region I) (Region H) 



a 

r“ 

*^Wi 


rji 


(Region I) 


(Region H) 


The net flux Fi from layer I at 
T = and from layer II are 


^ 2 .- - (M 3 + M 4 RS +M 5 Rw)' 


(AlO) 


t 5 


Vs = = -Tw^(M3+M4Rs+M5Rw)II 

(All) 


Sketch (b) Separation of shock 
layer (l) and injection 
layer (ll) . 

Dimensionless quantities in region II are based on wall (injection) tempera- 
ture T^, and boundary values and are to be determined for given 

Rg^ and Rg^^- With new definitions of 


Ri - Rv 


R2 = 


(A12) 


and energy balance at bomdary Tj_, 

^ ^2 = Ri ““ R2 


(A13) 


it follows that^ from eqmtions (AlO) , (All), and (AI3) , equivalent boundary 
values are 


I _ (M 3 + M4Rs) (1 - Ms^^) + T^^(Ma + M 4 Rs) 
Ri = R^ 


II 


- (Ms^ + Ms^^) 


(Al 4 ) 
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(Ms + M4 Rs)V V + (1 - M 5 I) (% + M4Rs)^^ 

- (Ms^ + Ms^^) 


(AI5) 


With these boundary conditions, flux toward the free stream, fg, and fliix 
toward the surface, from the whole system are calculated by 


fg = (Mq + MiRs + MeR^)^ (AI6) 

^wb = = -V(Mo + MiRs + MaRw)^^ (Al?) 


The radiative blockage fmction is then 


%b “ 


^ wb 


wo 


(91) 


where is the radiative function at the wall without injection (or 

ablation) and is written 


= (Ms + M4 Rs + MsRw)^ 


(AI8) 
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APPENDIX B 


THERMAL FUNCTION 


From thermodynamic charts or tables (refs. 3^ and 35)^ temperature varia- 
tion with enthalpy for constant pressure (or constant density) can be 
expressed around reference temperature T^^ as 


T^ = 




^r) + ^ 


(^) (h - 
/r 




(Bl) 


For the present work, only a few terms in eqmtion (Bl) are sufficie nt t o 
cover the range of temperatirre of interest; however, average slope (dl^/dh) is 
effective and simple for general applications. Thus^ with temperatiore behind 
the shocks Tg^ as reference^ 


t4 




hg) 


This equation, normalized by dividing by Tg"^ becomes 


(B2) 


1 - a(Fg - F) 


Where 



(B3) 

(Bi^-) 


This equation can be called the ideal or linear thermal relation since the 
transfer equation is linearized (in contrast to the linear relation between 
enthalpy and temperature for an ideal gas ) . 

The temperature relation for constant density in the preheating zone can 
be written as 


where Tf is the temperature just ahead of shock front and is to be deter- 
mined; Tf is selected as reference temperature rather than free-stream 
temperatiure T^ since transport properties depend strongly on temperature 
near Tf . In normalized form^ with Tf as reference temperature^ 
e qmt ion ( B5 ) bee ome s 


53 





(T/Tf)^ 1 _ a^(F^ - F) 


(b6) 


Temperature can be calculated from the relation derived for other refer- 
ence gas temperature, as follows: from equation (87b) and for h|- - h;^ 


where 


h™ - h 

r (x» 




Mq(1 - 2E3T^) 

(r - r ) = 

1 - oMo(l - 2 EsT^) 


(B7) 


(b8) 


Equation (b 8), which takes only absorption effects into account and thus 
neglects gas emission in the preheating layer, is calculated from equa- 
tions (83) and (86) for given and Too- The value of temperature Tr^ which 
corresponds to enthalpy h;^., can be obtained from the real gas chart for con- 
stant density (p* ^ p^) . Since Tf - T^,, one can use (8T^/8h)j, for any 
temperature near T Tj so that in the preheating zone 

T/ T/ + (hf - hr) (B9) 


After equation (B7) is substituted into (B9) and divided by T 

s 

T/ T/ + ^[(Bsp - r^)' - (fgp - r^)r] (BIO) 


Where Fgp and including the effects of both self-ab sorption and gas 
emission^ are now calculated for the same and as used in equa- 
tion (b8 ) . The thermal function is used to calculate the radiative flux with 
preheating and is presented in appendix G, 

With this temperature Tp and d^ ^ the temperat-ure relation given in 
equation (b 6) can be used to evaluate radiative transfer in the preheating 
zone by a procedure similar to that given for the shock-layer analysis. 

Similar thermal relations can be used for the injection layer except now 
constant pressure over the layer is assumed; thus^ 

(T/T^)^ 1 - d^(F^ - F) (Bll) 

where is defined by equation (99), and wall temperature, T-j^, is the 

reference temperature for normalizing the thermal function. 
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Note that the linear thermal relation may not apply in some gases; hut it 
can he used to’ predict the importance of emission effects (from the preheating 
and the vapor layers) on radiative heat at the wall. 


Temperature-Pressure Relation 

The temperature variation for other pressures (fig. 2) can he written 
empirically as 


= T^(Ps,h) + [1 - T^(Ps,h)][l - (P/Pg) (B12) 


where Pg is reference pressure near p^ and m is constant; for example^ 
m - 0.69 for Tg ~ 15^000*^ K and for pressure ranges of Pg = 10“^ to 
10 ^ atm. 
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APPENDIX C 


MODIFIED BESSEL FUNCTIONS OF THE FIRST AND SECOND KIND (l^K) 


The differential equation of radiative transfer for stagnation flow 
(7 = 1) is (eq. ( 3 ^)) 


d^r 2mndT^ df - O 

de^ g de ^ 


(Cl) 


which can be reduced to the modified Bessel equation. The general solution is 


r = 6 [cily(n0) + C2Ky(n6)] 


(C2) 


Definitions and properties for these functions are presented in many textbooks 
(refs. 24 and 25)* The detailed calculation scheme involved in the present 
application is given in this appendix. Bessel functions in ascending series 
are (refs. 24 and 25) 2 K 


Iv(ne) 




K!r(v + K. + 1 ) 


(C3) 


K=0 


where r(v+K+l) is the gamma function. 


K (ne) = g . . A [I_v(ne) - Iv(ne)] 
2 sin Ttv 


(C4) 


Properties of the p-unction 


Differentiation formula Let the flux function r^(0) and ^^^(0) he 
defined hy 


r^(0) s ciz^(e) + C2z_^(e) 

r^je) = ciz^^(e) - c 2 Z_(^^^(e) (v ± i 7 ^ v) 


(C5a) 


where 
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z^(0) = 0 ^I^(ne) 
z_v(0) = e%(ne) 

Since derivatives of Zy fmctions are 

A [Zv(a)] = n9Z^.j(e) 

^ [z_,(e)] - -nez_(^.^)(e) 

it follows that 

^ r^(0) = n0[ciz^_i(e) - C2z_(^_i)(e)] 

or 

^ r^(e) . n6i„^(9) 

Recurrence relation *- The recurrence formulas of and are 

2vl^(n0) = n0I^_^(n0) - n 0 I^^ 3 ^(n 0 ) 

2vKy(n0) = -n0Ky_i(n0) + n0Ky+j^(n0) 

With slight modification of the last two equations it follows that 

2vr^(0) = n0^r^_^(0) - nr ^+ 3 ^( 0 ) 

From equation (C7) and (C9) ^ one obtains 

e A r^(e) = 2 vr^(e) + nr„^(e) 


(C5b) 

(C5c) 

(C6a) 

(C6t>) 

(C7) 

(C8a) 

(C8b) 

(C9) 


(CIO) 



therefore with eqmtion (CIO) it becomes 


— [er^(e)] = (2v + i)r^(e) +nr^^^(e) 


(cii) 


Integrals involving r,XQ) functions.- The integration involving the Fy 
function, which appeared in the evaluation of boundary conditions (eqs. (25b) 
and (c)), is 


±ne 


Integrating by part^ one obtains 


±ne , , r ±ne , 

I e r^_^(0)d0 = 0e - / 

O L Jq J O 


V +ri0 

e [±n0r^,_i(0) + 0r^_^(0)]d0 


where 




From eqviations (C7) and (CIO) it follows that 

e*'"®[±n0r^_j0) +0r^_^(0)]d0 

^ ^ oT ±n0 

= 2(v - 1) / e r 


w 



±n0 ' 

fo ^ 1 

e r^(0) 

_ _ 


d0 


or 


w 


±n0 


e rv_i(e)d0 = — i e 


±n0 


2v - 1 


-l''^W p 

er^_3^(e) + r^(0)J Q 


Rv > i 


Finally, it becomes 


' ”^w ±n0 , . 

e r^_^(0)d0 = 


1 L±nr 


2v - 1 


w 


TwRv-i('^w) T Rs 




(C12a) 


58 



Then integration associated with Fg in equation (25t>) is 


f ® - 'r)dT = I Vv-I^'^w) - ^’s + e (C12b) 


and that of F in equation ( 25 c) is 


/ e - T)dT = |e TwF^_^(t^) + Fsj - F^| (C12c) 

Integral of the Z^(9) function.- By definition, Cy(©) is expressed as 


nd r9 

c^(0) s J Z^(0)d0 = J e^i^{nd)de 


Similarly^ 




0’^Kv(n0)d0 


The term C^_^(n0) in equations (37) are calciilated by 


V— X 

Cv-i(0) = I (l) r(|) r(^v - - Lv->0)iv-a(ne)] 


(C13a) 


also. 



l)® + L^_Jn0)K^_2(n0) 


(01313) 


where L^(n0) is the modified Struve fimction (ref. 25) for which a series 
expansion can he written as 


L^,(n 0 ) 


v+l oa 



(n0/2)^“^ 

r[m + ( 3 / 2 )]F[v + ( 3 / 2 ) + m] 


(Cl4a) 
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In the reciirrence formula 


L (ne) = (f ) ^ i — (l + y 

\ 2 y r(3/2)r[v + (3/2)] \ ' 


\ = 1 




(2m + l)(2v + 2m + 1) 
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APPENDIX D 


EVALUATION OF BOUNDARY CONDITIONS 


The eqimtions required for calculating the boundary conditions are 


^ - 2nmd (^ ) ^ - n^P = 0 

d02 \eJ d0 


and the boundary conditions are 


r. T, ^ T^ r ” si4 

Fg = -2inRg + 2mR^ + mn / T^e dt 


r 


= -2mRge + 2mR^ - mn T T'^e 


's'" 


The general solution for the flux function can be expressed as 


where 


r(e) = ciXy(n0) + C2Y^(n0) 
>CQ(n0) = and Yo(n0) = 


(7 = 0 


Xx(n0) = 0 Iy(n0) and Yi(n0) = 0^^(n0) (7 = 1 


The thermal function F is related to the flux function by 

dF 


dF = 


ie/r^) 


7 


and to the temperature by 


= 1 - a(Fg - F) 


( 25 a) 


( 25 b) 


( 25 c) 


( 26 a) 


) 


(16) 


(15c) 
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The boundary condition can now be evalmted, say for 7 = 1, from 
equations (25) as follows; 


r = mn 
s 


e (l - aFs)dt + mna J e F dt - 2mRg + 2mR^ ^ 


The second term on the right side of the eqmtion becomes 

r "'"w —at — nT,j r '*'w n0 , , 

/'^e Fdt = e ^ e F(0)d0 


= --e 


1 -nTv/ Jl0„^ ^ 1 ,-nTw r ^ ^n0 


e“®F’ (0)d0 


By eqmtions (36) and (C12) it follows that 


r = m(l - aF )(l - e + ma(Fg - e - 2mRg 


+ 2mR^e 


-aT-w 


mnar^ 
2v - 1 




+ e 



But from eqmtions (l5c) and 


V = mnaT^ + ^ (35^)) 

after the terms are rearranged it follows that 

“rx^ 

fs = m - me”’^’^'^[l - d(Fs - F^) ] - 2mRg + 2mR^~“'^w _ | [T^r^_^(T^) - fg +e 

(Dl) 

Similarly^ the bomdary conditions for with use of eqmtion (C12a) or 

(C12c)^ can be written as 

1 

= me"''"''^-m[l-a(Fg-F^)]-2mRge~’^’'''+2mR^ - | {e'^'^^T^v-iC'^w) + ^’sJ * ^w} 

(D2) 


Since equations (37), (C7), and (39) provide that 
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Tg - ciZy + caZ_^ , 


* 


- C2Z_(^_^) 

Fs — Fw ~ n''''w[ciCy_i — C2C_^y_j_)] J 


(D3) 


two simultaneous equations^ with Ci and as unknowns^ are obtained from 
equations (Dl) and (D 2 ) as written in equations ( 32 d) and (38a) . Similar 
procedures were applied for evali;iating boundary conditions for 7 = 0^ and 
these results were given in equations (32c) and (32d) . 
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APPENDIX E 


ASCENDING EXPANSION OF Z FUNCTION (t^ « l) 


For 7=1, ■when v is fixed and the argument approaches zero, the 
follo'wing ascending series of the modified Bessel functions are useful for 
evaluating fluxes for « 1. From equations (C3) and (C4) 


Zv(e) = e^iv(ne) - 


V 


r(v + 1) 


(l/2)ne 1 / / \ 

1 + — ... (v -1, -2, . . .) 

V + 1 J 

(Ela) 


=e\^(ne) 


2 sin jtv (_ 


- Zjne) 


- 1 (if V - • ■ •} 


(Elb) 


z!v ^ Z-v(O) = I (l) '■(•') 




^V-l g2(v-l) 

r 1 

' r(v) 

_2v - 1 [2/ + 1) 





- 2 (a-v) 


+ |H 


1 v(2v + 1) 


n 


03 


2/ 3(2 - v) 


(E2) 

(E3a) 


(E3b) 


With these supplementary forms, the qiiantities in equations (38b) are 
expressed by the notations: 


Qi - 



T 


- 2 (x-v) 

W 


(E4a) 
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Replacing 0 ly in equations (Ela) and (E3a ) , one obtains 


i-v 

% - (r ^ 


r(v) 12 V 


" T "w"" 0(V) 


nT■^J^(l - 2]maT )Qi 


i-v 

2 A 1 


^v-i Vn 


r(v) L 


-a(i-v) 


• w 


+ l“) T. + 0 (t/) 


27 V w 


1+1 


'v-x - Vn 


^ i-v 
2 \ 1 


r(v) L 2 v - 1 


1 2 V -1 / ON 

T.. + 0(t 2) 


w 


+ o(+/) 

. 2 mncL 


Q] 


(E^Td) 


(E5) 


(e6) 


(E7) 


From eqviatlon (Elb) one can write 

"-V -- [l (l)'’'-*''*]!" * - [i icha: (l)” •F(7^K'’‘" " 

by using the following properties of the gamma function (refs. 2h and 25 ) 


(e8) 


r(v + 1) = vf(v) , = sin_rtv 

r(v) « 

and 

l(i - ^o) ~ n/«(1 ^ ^0^0^ ^ 

where = mnar^ and cp^ = 2 Zn 2+7 (Euler’s constant). Consider the first 
bracketed term in eqixation (e 8 ); one obtains 
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1(1) r(v) 



2 V -1 


r(v) 


r(i - v)J 

(1 - 29oVo)Q2 


iQg 


(E9) 


Consider the leading part of the second term in equation (e 8 ) . One obtains 

I (ij ^ " (H) 

^ (1 - 2vo)nQ2 (ElO). 

Finally, by combining equations (e 8 ), (E 9 ) , and (ElO), one obtains 


Z_v - [1 - (2mnacp^ + nT^^^°)r^]Q^ 


(Ell) 


i-v 


Similarly, -^t he function follows 

7 ~ 1 ^ r(v) _^l/2\ ' V 2 (v-i) 

Z_(v_i) - - 2 vs) ■" 2 W 

- (1 - 


(E12) 


where 


1 

2 



r(v) 

1 - V 


^ n[l - 2(9 q - 1 )Vq]Q 2 


From equation (E 3 b), one could directly obtain 


^ ~ r _ i _ _ r,T ' 

'-(v-l) [ Pinn a ^ w 


2 ( 1 - Vo)' 


-2(l-v) 


w 




(E 13 ) 


Through equations (E5) to (EI 3 ) , K and M are reduced as given in 
equations ( 55 ) and ( 56 ) after several simplifications. 
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APPENDIX F 


ASYMPTOTIC EXPANSION OF THE Z FUNCTION (t^ » l) 


Uniform asymptotic expansions of the modified Bessel function for large 
argijments (refs. 2k and 25 ) are 


Iv(v0) - — — 




n/2«v (1 + 02) 


1/4 


1 + 



(Fla) 


KJ v0) if Tl - o(i^ 

y2v(i+ 03)1/4 L \vjj 


where 


n = -v/T 


+ + In 


1 + \fl + 0^ 


(Fib) 


(Flc) 


Evaluation of and K^_^(nT^) 


Since v is 


V = mnar^ + | 


nT„ ^ (v - k) 

»» mrw 


ma 


1 ^ k - . (V .?,). 

V - k 


(F2a) 


let 


z = and z 

^ mx 


. ji + w ?) : 

L V - k - 


(F2b) 


then 


7TT? « A + [i .(ic - 


(F2c) 


Similarly, 


1 + ^/l + 2^ - (l + 7 1 + 


^ ^ k - (1/2) yi + Zo^[Zo^/(l + Zq 2)] 
V - k 


1 +^^ 


+ Zr 
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and 


In 


s 


s: In 


zo 


1 + + Z" 


1 




V - k 


1 +Jl+ Zq^ n/i + zo^ 


From equation (Flc) it follows that 

(v - k)Tj - (v - k)/^ + Zq^ + In 


1 +^/^ 


+ z. 


(F2d) 


+ ( k - i Y 1 - ^ , . _ 

Yl + Zq 1 + + ^o^ 


+ z. 


- (v - k)qQ + (^k - 0 Ho 


(F3a) 


where 


Therefore 


%=J- 

Ho =«yi + Zo' 


1 + Zq + Zn 


Zo 


1 +JT7^ 


= \Xq - In Pi 


[k-(i/a)]Ho , , 


(F3h) 


(f4) 


How, Pi and Pa from equations (28) can he also expressed as 


(F5) 


Using equations (f 4) , (Fla), and the definition in (F3h) , one can obtain the 
following asymptotic recurrence formula; 


Iv_k(“"^w) - Pi 


(f6) 


Similarly, for as one obtains 

-[k-(i/a)]Ho„ / 

Kv-k(^-^w) - ® 


(F7) 
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The recurrence formula thus becomes 




(f8) 


Asymptotic Expansion of 


As given in equation (C13a) 
v~i 


= i (l) Xf) - i) ® Iv-i(“0)Lv-2(a0) - Ly_3^(n0)l^_2(n0) 


The modified Struve function L^. (ref. 25) for 0 » 1 is 


Lv(n0) - l^y(ne) V — 

Z_i p[y ^ 


«l 

(-1) r[k + (1/2)] 


v 2 k-v+l 


r[v + (1/2) - k](n0/2)‘ 

(|arg 0| < n/2) 


(F9) 


(FlOa) 


Let 


’'k,v 


k+1 

(-1) r[k + (1/2)] 


r[v + (1/2) - k](n0/2) 


2k~V+l 


then 


^k,v 


-[k-(l/2)] [v + (1/2) - k] 
{nd/2f 


7 


k-1, V 


(FlOb) 


As V and 6 v ^k-^i y ’ "therefore, if is siifficient to take only the 

leading term in equation (FLOa; as 


Ly(n0) - I_y(n0) 


1 (n0/2)^ ^ 

r(i/ 2 ) r[v + ( 1 / 2 )] 


(FUa) 


With equation (FUa) the bracketed qmntities in equation (F9) can be 
expressed (for simplicity, argument n0 is omitted) as 




-i^v-a 








V-1 


I 


V-2 


^-(v-i)I 
“ 7 o,V-2^ 


v-aj 

V-J 


(FUb) 
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Since Wronskians (ref. 25) vanishes as 0 -> <» 


V-2)J [^V-l^-(v-2) ^-(v-l)^V-2j 


— ~ sin rt(v - 2) « 1 
rt n0 


It follows that 


I V- V-2 ^ V- V-2 J 


1 (n0/2)^~^ 

r(i/2) r.[v - (1/2)] 


[v- (3/2)1^ 
>-2 'L ne/2 


If d = and 00 , then, from equation (f6) , 


■ _ V - (3/2) 1 . 

1 '"“ nT„/2 >'-^J •' 

Finally, equation (F9) , with equations (Flld) and (File), becomes 

Cv_i(tw) ~ ^ V Iv(’^'^w) = 


From equation (E3a) the same resiilt is obtained for « 1 and v » 1 since 


Cv-l('Tw) ^ 


Zy 


- (l/2) riT^ 


An evaluation similar to the above procedure can be applied for (©) 

except that one substitutes Instead of I^j the results are 

^|^^v-i-’^-( v-i) “ ^-(v-i)^v-2 ^-(v-2)^v-i ~ n0 


1 (n0/2)^~^ 

r(i/2) r[v - (1/2)1 


Kv_2 + 


(3/2) 




(0 oo) 


which lead to 
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This approximation can also be compared with the exact form of = oo in 
integral tables (ref. 25 ). 


J thCy(t)dt = 2^ ± v) >-i 


(F13b) 


Note that substituting p = v in eq.mtion (F13b) gives the approximate 
formula fomd in equation (Fljar) . Note also that if v » 1, but « 1, 
the following result is obtained after integrating equation (Elb) 




nT. 


w 


2(v - 1) 


-V 


(Fl^l-) 


Flirther asymptotic expansion of the following functions are necessary 
to simplify the asymptotic evaluation of the flux function. Thus, as v oo, 
the gamma function can he expressed as 


and 


, X y — -V v-k-(i/2) 

r( V - k) - e V 


lim I 1 +^ \ -* 

V>oo 


(F15a) 

(F15b) 


It follows that 


= I (i) 


/ire /^2mg 

72 v\e w 


(F16) 


Equation (F13) then becomes 


From equation (Fla) 


^-(v-i) 


^ (l/2)r(l/2) 


-i^oA 


-V 




/ ’’lo \ V 

(e T^) 


(FIT) 


^ Oae"%)" 


j2nvy.^e^o 


(F18) 
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Also from equation (f 6) it is shown that 






' w 


From equation (Flh) one obtains 

^-V - 


and from equation (f8) , 


^ Pa^-v 


With the asymptotic values derived above, the following relations hold; 


Since 


mnax^^_^ 


^ i = i (1 + Pi)(l - Pa) 
zo 2 


mnaTwC_(^_3^) 



-nT-^ 

e 


-[v-(V 2 )]/nia 2 io /2 --Zq V 

e e {e 


one obtains 





l+Zf 


Zq® 




Then 



since - 1 > 0. Also one finds 


-nr^ - 

e mnaT^_(^_ 3 ^) 


Z-v 



I+Zo-Po 

V 


2/(1 + IIq) 


l+Zo-lio 


« 1 


(F19) 

(F20) 

(F21) 

(F22) 

(F23) 

(F21)-) 

(F25) 

(F26) 

(F27) 
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since 1 + Zq - [Iq > 0 . And 


- 




^-v 


^^0-1 VI + 


Mo-i 


Uo 


» 1 


(F28) 


From equation (F28) it is obvious now that 


mncc/T ^ ^ ^ ^ 



» 1 


(F29) 


The values Kq and (eqs. ( 6 l) and (62)) are readily evalizated from the 
asymptotic values presented above . 


Calculation of Radiative Flux Function f 

Coefficients ci and C 3 of the radiative flux function^ f, can be 
evaluated from the previous results; thus^ 


Cl 


-kgke + k4ks 
A 


Ml - P2)[(l/2) - Rw]Z_v + MV^ + l)[(l/2) - 

(1 + Pi)(n^ + 1)Z^Z*^ 


4m[(l/2) - Rg] 

(1 + Pi)Zv 


(F 30 a) 


Similarly, constant Ca is evaluated as 
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kilss - ksks 

C2 = 

A 


^ -Ml + Pi) [(1/2) - R „3z^ + (8/zo)m[(l/2) - RgJZ^ 

(1 + 3i)(N/lr? + 1)Z^% 


-2m[(l + P2)/(l + Pi) - 2(1 - p2)Rs + 2RJ 
(<77^ + 1)Z*^ 


(F30b) 


Since the flux function is given as 

r^(0) = ciz^(e) 

it is necessary to examine the asymptotic 


+ C2Z_^(e) 

behavior of the functions 


Zv 


and 


z_v(e) 

zfv 


(F31) 


for 1 « T « Tw Similar evaliiation made for eq.\jatlon (F^J-) may be applied 
for obtaining the asymptotic form of I^(n0). Thus, 


n0 = nr^(l - t) = vZqCi - r) (i - ^ 


that 


Iv(iie) “ e 


-(l/2)H* 


I^(vz*) 


that for equations (F2c) and (F2d) . 

[i* = Jl + (z*)^ u ^1 — - t') 

°\ 1 + Zp2 y 

7 ^ 


• — ) 
2vJ 

(F32) 

(with 

k = O) by z* so 


(F33a) 

manner 

analogous to 


(F33b) 



Zn 


(F33c) 




1 + M* 


(l . r.lng^-^ 

\ ^^o / ^ 1^0 


then 


n* - (^o “ Pi) 


- (fof + i) 




Thus , from equation (Fla) it is shoTO for 1 « t « that 


Iv(ne) (l +i 


1 o 


2 \ -(l/^)(l/^o)T “liQinnaT 


2 1 + Zo^ 


T e 




-HonmaT 

- e ly(nTw) 


2iy(0) _ V _|j,Qii]n,a,T -(i-Hio)nmaT 

“(l-T)e =e 


-nPiT 


- e 


Similarly, Ky(n0) can be written as 


li^mnaT 

K (n0) e K^(nT^) 


and so 


Z_v(0) , 

— =r (1 _ t) e 


^ ^ ^ ~ e 




Therefore, from equation (F26), it becomes 


Z-v(®) Z_y(0)/z_y -1/2/1 + “UPaC^’w"'’’) 


* -,* /r, ~ *^0 V 2 

V 


zrv z*vA 


9 


(F33d) 


(F3i^) 


(F35) 


(F 36 ) 


(F37) 


(F38) 
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For sufficiently large Ty the value of equation (F38) is negligible compared 
to that of equation (F35)> for example, 


«l 



(F39) 


Thus^ the asymptotic radiant flijx fimction becomes 


r^(e) « ciz^(e) 


2m 

1 + Pi 


(1 - 2Rs)e 


-nPiT 


(f 4 o) 
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APPENDIX G 


CALCULATION OF RADIATIVE FLUX WITH PREHEATING INTERACTION 


The energy equation in the preheating zone is^ from equation (77)^ 
Pco^co^S A* - 


2oTs^V 




(Gl) 


It follows that 


h. 


- - rj 


(G2) 


vhere 


r* = q.*/2rtB 


(8a) 


Defining Fgp as the radiative flux of the preheating zone at the shock vave 
and hf as enthalpy just ahead of the shock fronts one gets the relation 

\ =" Aa(Fsp - r^) (87b) 

S 

Preheating increases the temperature behind the shock by the amount 

Tgp^ = 1 + a(F3p - r^) (80) 

Flux (coming from the hot shock layer) at the shock front is 

^sp = (G3) 

where R 3 '^ = and R^"^ = R^/T 3 p^- With equation (8o) it follows that 

Tsp = [1 + a(F3p - rj]% + M1R3 + MaR^ (G 4 ) 

The continuity of the flux function requires that (see sketch (a), main text) 

I'sp = % - »s (05> 

The temperature ahead of the shock T^ is determined by 
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(bio) 


'^r ” ^00 ~ (^sp ” 

where is the reference temperature derived from the assTomption that only 

absorption by media is a major factor in the preheating zone (see appendix B, 
eq.s. (b 8) and (BIO)). Thus, boundary fluxes of the preheating zone from 
equations (26) are 


= ff^[Mo + MiRf + Msi^]* 

(g6) 

= Tf^[M3 + M4Rf + Mgl^]* 

(G 7 ) 


where superscript * designates quantities in the preheating zone, and 

Rf* = 


It follows that 


^sp - ^ 0 ° = -Tf^(Mo - M 3 )* - (Ml - M4)*Rf - (M 2 - Mg)*R^ (g 8 ) 


Now, one can solve the simultaneous equations given in equations (G 3 ) 
through (g 8) for Boo = 0, R^ « 0 and for given optical depths (too and t^) of 
the preheating zone and shock layer, respectively, as 


-T ^ 
g 


Mq *(1 + M *) - Ma*(l + Ml*) - 


(Mq- M3) * 

1 + Ml 


sp 


Too = 




+ (Ml - M4)* 


Mq 

1 + Ml 


A 


(G 9 ) 

where 

V = - 2 E 3 (tJ ]/[1 - oMo[l - 2 E 3 (too)]} (GIO) 


and 


A = l+Mi*- 


1 _ 


1+Mi 


' CC/j^ 


Mo*(l+ M4*) - M3*(1+ Ml*) - 


Mo*-M3*- 

1+ Ml 


-a(Mi- M4)* 


Mq 

1 + Ml 


(Gll) 


Then Tgp'^ and T^^^ will be obtained from equations (80) and (BIO), respec- 
tively. Thus, T^gp'^ may be written in terms of known quantities as 
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1 r 

(l + Mi*) - + |^Mo*(l+ M 4 *) - M3*(1+ Ml*) - 


T ■* = 
■^sp 


(G 12 ) 


Boundary -values are given "by 


Rf = - 


^ - ^3*) 

Ml* - M4* 


(GI3) 


Rs = 


Rf 


[1 + g(rsp - rJ]Mp 

1 + Ml 


(G 14 ) 


The terms Tgp and can be obtained from equations (g 6) and (G7) or are 
simply 


Tgp = Rf - Rs (GI5) 

^co = Rsp - (^sp - r^) (G16) 

From equations (GI5) and (G 16 ), the radiative leakage parameter, vhich now 
includes emission effect, hecomes 


^00 ” ^ooAsp 


(GI7) 


Note that if the temperature jump due to preheating is not negligible 
(Pgp “ Pg 7^ 0), the following modifications are made to calculate the enthalpy 
and temperature distribution through the shock layer* Quantities K5 and Ke 
in equation (38a) are replaced by 


K5 = 2m 


Ke = -2m 


I (1 + 5Arg) 

^ (1 + ccAFs) 


- RgJ - 2me'‘^'^w ^ (i + 

- R^j + 2me"“'^w [l 



(G 18 ) 

(GI9) 


where APg = Pgp - P^, However, the same result will be obtained, as one can 
see directly from the last equation, if the boundary values Rg and R^ are 
divided by Tgp^. Thus, from equation (G 3 )> we define. 



s 


1 + oAP, 


s 


(G 20 ) 
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and the thermal function (Fgp - F)"^ as given hy equation (88a). The enthalpy 
and temperature distribution given as eq.uations (89) and (90) will yield 
results identical to those derived from equations (GI8) and (G19) since all 
other K in equation (38) are constants for given a (constant) and 



APPEiroiX H 


CALCULATION OF RADIATIVE FLUX IN INJECTION AND ABLATION LAYERS 


The general solution of radiative heat transfer for the two -layer prohlem, 
shock and injection layers^ for example, has heen presented in appendix A. 

The details of the analysis for both absorption and emission of injected gases 
on radiative heat to the wall is considered here. Since the characteristics 
of absorption coefficients for injected or ablated gases are not yet well 
known, the absorption coefficients simply appear in this analysis as a 
parameter, vhich is independent of temperature. 


Injection Layer 


For radiative transfer with gas injection from the body surface, injec- 
tion rate, is an arbitrary parameter and is thus independent of the heat- 
ing rate at the wall. For this case, the analysis of radiative heat transfer 
is straightforward. 


Let flight conditions be fixed and assume the optical thickness 
(or in the shock layer to remain unchanged during injection since the 

shock-layer structure is relatively insensitive to the presence of an injec- 
tion layer. Note that assuming a given value of is equivalent to assign- 

ing a given body radius (see eqs. (ylib) and (76a)). The following outline 
illustrates the calculation procedure: 


(a) Assign (t^ and a) in both regions: Region I (shock layer); for a 

given body radius (and flight condition) one can find and a. Region II 
(injection layer); optical thickness, T^'b, is assumed. The injection rate, 
f-vT-, is also assumed, and one finds 05 from equation (99). 

T TTr 

(b) Calculate and from equations (a 14 ) and (AI5) . 

(c) Calculate the absorption coefficient, which corresponds to 

quantities (x^, a, Rg, by equations (7^b) and (95b): 


(pOs 




(pK)sks 


J\ 


Ps 


M 

00 

Mb 


f-j^R 


(HI) 


where 

(p'")b^o = '^yrb f (P/Pw)""^^ (Ti*) 

’-'o 

(d) Calculate and from equations (AI7) and (9I) • 
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II 


II III I III II I II I null I III I 1 1111111 II II nil nil I III II null I nil I HI l■H■ ■ ■■■ mm 


By repeating procedures (a) through (d) for other values of and 

ve can obtain radiative hlocbage function as a function of the injection 

rate with the absorption coefficient ratio (pK)-| 3 /(pK)g as the parameter for a 
given hody radius. 

It is interesting to note that has nearly a single correlation curve 

as presented in figure 10(a) . 


Ablation Layer 


For an ablating surface^ the ablation rate^ f^, is coupled with radiative 
heat at the wall Thus, some simple algebraic iteration is required to 

solve equation (98)» Using the approximation given by equation (92), one can 
calculate on initial input value of from equations ( 98 ) and (99)# by 


% - 




^ ^ VO 


2E3(t^) 


(H2) 


for given Similar procedures (as outlined in the previous section) are 

then used to calculate and r-v/j-,. This provides a new abla- 
tion rate, and thus, previous calculations are repeated until the value of 

f^ converges. Only two or three iterations are necessary for good conver- 
gence. Even the results obtained by using equation (H2) as an initial input 
show siifflcient accuracy. With this ablation rate one can compute the absorp- 
tion coefficient ratio from equation (Hi) . Thermodynamic properties of carbon 
are obtained from reference 36. 
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APPEiroiX I 


EQUILIBRIUM AIR PROPERTIES 


Thermo dynamic charts are presented for flight speeds to 30 km/sec in 
figure 13 . These charts give the pressure and temperature behind and the 
density ratio across the normal shock wave, and can be used to obtain refer- 
ence properties for the present solutions. These properties were computed 
from the thermodynamic charts of references 3^ anti 35* 

Radiative properties for equilibrium air are presented in figure 14. 
Figure 14(a) presents E^/2, one-half the emission rate per unit volume, for 
optically thin air. This chart was prepared by corrbining the calculations of 
reference 37 for low temperature and reference 38 ^*or higher temperatures. 
Figure l4(b) presents Planck mean absorption coefficients from the same refer- 
ences used in the present analysis. The relation between emission rate and 
absorption coefficients is E.^ = ^Ps*^s^s'^* absorption coefficients 

for equilibrium air are not well established over the whole temperature range. 
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TABLE I.- EXPONENTIAL CONSTANTS, n 



n 


n 

0 

2.000 

.40 

1.736 

.002 

1.990 

.50 

1.709 

.004 

1.986 

.60 

1.686 

.006 

1.982 

.70 

1.667 

.008 

1.978 

.80 

1.651 

.010 

1.97^ 

.90 

1.636 

.012 

1.970 

1.00 

1.623 

.014 

1.966 

2.00 


.016 

1.963 

3.00 

1.519 

.018 

1.959 

4.00 

1.507 

.020 

1.956 

5.00 

1.503 

.040 

1.929 

6.00 

' 1.501 

.060 

1.907 

7.00 

1.500 

.080 

1.889 

8.00 

1.500 

.10 

1.873 

9.00 

1.500 

.20 

1.812 

10.00 

1.500 

•30 

1.769 

00 

1.500 
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Figiare 3 *- Comparison of exact and approximate solutions 
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(b) Enthalpy distribution in stagnation flow. 
Figure 3.- Continued. 
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Figure 7*“ Distributions of temperature, radiative heat flux, and heat-flux 

derivative in shock layer. 
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Figure 9.- The effect of nonadiabatic flow on shock-standoff distance. 
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